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Sicut aquae tremulum labris ubi lumen ahenis 
Sole repercussum, aut radiantis imagine Lunae 
Omnia pervolitat late loca jamque sub auras 
Erigitur, summique ferit laquearia tecti. 

[As when ruffled water in a bronze pot reffects the light 

|Jy»|2i"f 5i C*CX of the sun and the shinning face of the moon, sending 

V.^ J.C>i.V> v.^ shimmers flying high into the air and striking, against 

the paneled ceilings. Virgil, Aeneid, VIII, 22. ]El 



We are not busy with some particular physical prob- 
lem. We are "ranging to and fro over the wastelands 
of" the modern physics similar to Virgil's sunlight reflec- 
tions. This thesis consists from a number of parts. They 
are completed separate works motivated by recent ex- 
periments. The words "superlattice" and "current" are 
common in these works and widely used. Two chapters, 
1 and 3, were written especially for this thesis. They are 
actually unpublished results and discussions. The kind 
atmosphere of the Racah Institute of Physics helps us a 
lot. I'm grateful to this place. I would like to acknowl- 
edge my and Boris Laikhtman's discussions with Hol- 
ger Grahn, S. Luryi, Y. Lyanda-Geller, Yehuda Naveh, 
M. Raikh, Leonid Shvartzman, and V. Zevin. I thank 
F. Bass for permission to use Part II, "Quantum super- 
lattices" , of his book, prior to publication. 

All the principal results of this work concern the verti- 
cal transport in generic three dimensional superlattices. 
In the 1st chapter we make an historical introduction, 
then we discuss the geometry of the problem, and the 
physical parameters associated with structure of electron 
minibands and strength of external fields. We found the 
effect of collisionless transverse magnetoresistance, and 
we discuss it in the 2nd chapter. This effect is similar to 
collisionless Landau damping in a plasma and we utilize 
the same name. In the 3rd chapter we provide quantum 
mechanical reasons for the above effect; we show how a 
magnetic field bends narrow superlattice minibands, and 
we classify the states into Landau-type and Stark-type. 
In the 4th chapter we compute longitudinal magnetore- 
sistance of superlattices due to the imperfections of the 
interfaces. Correlation length of the interface roughness 
can be measured independently by this method. In the 
5th chapter we discuss the current-voltage characteristic 
of superlattice when an electric field destroys the one- 
miniband transport. We found that the structure of the 
high-field domains in the superlattices is complicated, 
but can be described analytically with great accuracy. 
This structure reveals itself in the details of the current- 
voltage characteristics. All our results are consistent with 
existing experiments, and we make careful comparison 
of theoretical predictions and experimental results in all 
chapters. 
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Chapter 1 

Introduction 



1.1. GENERALITIES 

We consider theoretically semiconductor superlattices. 
A superlattice is a one-dimensional periodic potential 
in monocrystalline semiconductor. The superlattice po- 
tential would be obtained by a periodic variation of al- 
loy composition introduced during epitaxial growth. It 
would consist of alternating ultrathin layers of two semi- 
conductors that closely match in lattice constant. 

We put forth a more specific definition: A superlat- 
tice is a stack of coupled quantum wells. Each quantum 
well limits electron motion in one direction and leads to 
quantization of energy; that is the elementary example 
in quantum mechanics textbooks. Coupling of quantum 
wells means that the electron can tunnel through the bar- 
rier separating adjacent quantum wells. 

Our definition is more accurate because variation of 
alloy composition leads to variation of all band structure 
parameters. They cannot be described by a simple peri- 
odic potential. They modify electron wave functions in 
potential wells and tunneling amplitude through poten- 
tial barriers. Our definition is particularly useful when 
variation of superlattice periodic potential is much larger 
than mean kinetic energy of electrons. In this case elec- 
trons are really sitting in the wells of the superlattice 
potential. 

The evolution of molecular beam epitaxy has allowed 
access to superlattices. High-quality devices started to 
be available to experimentalists in the early seventies. 
For example, effects of energy quantization and iunneling 
were observed in 1974, see review of Weisbucho 

Original motivations of superlattice research were to 
obtain devices with negativej-jdifferential conductance^, 
and amplification of radiationQiJ. These effects were pre- 
dicted for a superlattice placed in strong electric field 
perpendicular to the layers (vertical transport). Unfor- 
tunately, instability of electric potential profile in vertical 
(growth) direction and formation of high- and lowJield 
domains made realization of these ideas impossible.B 

Structure and dynamics of the high-field domains in 
superlattices is a separate lyanch of research. Attempts 
to buildfja microwave laserQ and a microwave radiation 
detectorQ motivated the investigation of electric field do- 
mains. The first project was not realized for many years 



due to the relatively small life time of excited carrier^. 
Progress in making a quantum cascade laser was reported 
only recently.E2l Structure of the high-field domains, i.e. 
potential profile in a vertical direction, remains to be a 
subject of intensive research. This potential profile, for 
example, carries Mot of information concerning the qual- 
ity of the deviccE^ 

There are some similarities between instabilities of the. 
potential profile in superlattices and in Gunn diodes.E3 
Therefore, one hopes to see high-field domains running 
through the superlattice. This was one of the possible 
explanations of the time dependent oscil.laijinQa.of the po- 
tential profile which had been reported.E3o'E3 However, 
numeric simulations show that the strengths of the high- 
and low- field domains change periodically in time with 
the domain boundary being pinned within a few quantum 
wells.ye 

Modern technology uses superlattices for construction 
of optical devices. The physics of such devices is tightly 
related to properties of excitons in quantum wells. We 
will not touch the world of optics in the present work. 
Besides the device application, growth of superlattices is 
still a challenge for labrotaries employing molecular beam 
epitaxy. 

Typical superlattice geometry is shown in Fig. 1.1(a). 
The horizontal sheets schematically show a superlattice 
sandwiched between materials providing electrical con- 
tacts. The superlattice itself is formed by alternating 
layers of wide and narrow band gap materials. The real 
potential profile of such a device is complicated, and de- 
pends on the materials and doping level. The part of 
the potential profile related to the superlattice is shown 
schematically in Fig. 1.1(b). The potential energy of elec- 



trons changes periodically in the y direction and remains 
constant in the x and z directions. 

This potential defines the sequence of the quantum 
wells separated by the barriers. Potential wells limit 
electron motion in the y direction and lead to energy 
quantization. We have chosen the energy of lowest elec- 
tron state to be the zero of energy. We show the case of 
two electron states in a well, they are marked by black 
rectangles in Fig. 1.1(b). The energy gap between these 



two states is Eg, and it is one of three main parameters 
characterizing the superlattice potential. 
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FIG. 1.1. (a) - Measurement of the vertical current in a 
superlattice. Current goes in the y-direction, transverse mag- 
netic field is applied in the z-direction. (b) - One miniband 
transport through a superlattice. The Fermi energy (dashed 
lines), Ef, lies slightly above the lowest level in the quantum 
well. Energy of the lowest level is zero. Energy of the second 
level is _E„ 3> -Ef- 



The superlattice structure brings to semiconductor 
physics a new length scale, which is the period of the 
superlattice, I, and a new energy scale, which is tun- 
neling energy, A. These are two other parameters men- 
tioned above, and their typical values in 1974 were 500A 
and ImeV correspondingly. Modern technology allows 
growth of superlattices having very fine structure. For ex- 
ample, the barriers of only three monatomic layers width 
were reported in Ref. |l^. The typical superlattice po- 
tential period in 1996 is about 50A and tunneling energy 
may reach 40meV. 

The resonance between energy levels in adjacent quan- 
tum wells leads to formation of minibands, and one can 
associate the tunneling energy with a quarter miniband 
width. This-energy can be measured by photocurrent 
spectroscopyllj and results arejii agreement with predic- 
tion of Kronig-Penney model E2l However, calculation of 
minibands of holes is a much more complicated prob- 
lem, beftaiuse wave functions of holes are four component 
spinors.u'O We will not consider holes in this work. 

The energy gap Eg is assumed to be very large in the 
theory of the one- miniband transport, much larger than 
the Fermi energy, temperature, potential drop of electric 
field per period, energy unc erta inty due to scattering and 
miniband width. In Fig. 1.1(b) we showed the Fermi 



energy between two minibands, i.e., the second miniband 
is empty. A high electric field destroys these minibands 
and leads to formation of low- and high-field domains. 
Inside the high-field domain Eg is of the order of the 
electric potential drop per period. For this reason, the 
physics of high-field domains stands separately from the 
physics of one- miniband transport. 

The tunneling energy is usually so small, that the co- 
herence of tunneling through subsequent barriers can be 
easily destroyed by external fields, sample imperfections 
or impurities. If it is known that sequential tunneling 
is incoherent, one can introduce the conductivity of the 
barrier jOj the transition time of an electron through the 
barrier We will use this approach in Chapters Hand||. 
If it is known that the tunneling is coherent, one can write 
down the usual kinetic equation for an electron gas with 
anisotropic dispersion. This dispersion contains all in- 
formation about quantum mechanical effects(tunneling), 
and the kinetic equation is already a classical object. We 
will use it in Chapter ^ Very often, it is not known 
whether tunneling is coherent or not, and it is possible 
to scan all intermediate situations by changing external 
fields. For this reason we will consider the quantum ki- 
netic equation derived by Keldysh's technique. It is the 
only tool which can describe the destruction of a mini- 
band by a high electric field. We will derive it in Chap- 
ter |. 
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FIG. 1.2. Measuredtfl I-V characteristics of the superlat- 
tice. The first discontinuity is the instability point of the 
one-miniband transport and corresponds to the formation of 
the high-field domain. Further increasing of the bias leads to 
the expansion of domains and small current oscillations near 
the value 7*. This is the first plato, see details in Chap. 

1.2. MOTIVATIONS AND THE PROGRAM 

This project is actually a collection of a few works 
concerning vertical transport in superlattices which were 
motivated by recent experiments. Motivation for each of 
these individual works and related experiments are con- 
sidered in their introductions, but we need to put all of 
them in the proper context here. 

A lot of experiments and theoretical problems concern- 
ing one-miniband transport were discussed in the liter- 
ature in the years 1970-1985, however interpretation of 
the recent experiments with high-quality devices is not al- 
ways possible in the framework of the old theories. Most 
of the recent experimental results are collected in the 
book of GrahnEa. The modern theories of superlattice 
physics are considered in this book only briefly. 

The simplest transport measurement is detection of the 
I-V curve. The typically observed I-V curve is very non- 
linea r, a nd it exhibits a sequence of discontinuities, see 
Fig. 



1.2 



The first discontinuity corresponds to the de- 
struction of single miniband. The nature of this instabi!-, 
ity point was considered for_Lhe first time by LaikhtmanEI 
and Laikhtman and MillerEj. The latter work contains 
the derivation of the quantum kinetic equation which 
is valid near the first instability point. The instability 
development leads to high-field domain formation. At 
least two states in each quantum well become involved in 






FIG. 1.3. Superlattice potential is shown schematically for 
several values of applied bias: before the first instability point 
(a), at the formation of the high-field domain (b), correspond- 
ing to expansion of the domain (c). 



transport, see Fig. 1.3 
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The structure of the high-field domain was calculated 
for the first time by Miller and LaikhtmanEZl. We repro- 
duce this work in Chap. ^. Since this publication, the 
structure of the domains has been intensively considered 
in the literature. For example, additional effects, like im- 
perfections and fiuctuation ot-barrier size, were consid- 
ered recently by Schwarz et alllil. Chapte r c onsiders the 
parameters of the typical 1-M.xupie, Fig. L2, in great de- 
tail. Recent measurementscaO'O are also in agreement 
with this theory. The physical problems considered in 
Chapter || are 

• limiting of the current by diffusion, 

• minimal domain size, 

• nonresonant tunneling of electrons inside the do- 
main, 

• time resolved process of the domain formation near 
the instability threshold, 

• mechanism of the domain expansion, 

• possibility of domain formation near the cathode 
and therefore injection of electrons into the second 
miniband and two-miniband transport. 

It is very well known that transverse magnetic field 
shifts discontinuities of the superlattice I-V curve to 
higher bias. Concerning the first instability point sue can 
mention the recent experiment of Aristone et alu ajod 
much more data can be found in the book of GrahnE^. 
It is believed that this discontinuity corresponds to the 
maximum of the I-V curve of the uniform superlat- 
tice. The position of the current peak is usually cal- 
culated by numeric solution of the kinetic equation in 
one mini hapd Hj Analytical results were obtained first by 
EpshteinEill23. They show a quadratic shift with mag- 
netic field of the I-V curve peak. However, when mag- 
netic field becomes stronger, the peak shifts linearly with 
magnetic field. The explanation was given by Miller and 
LaikhtmancJ in terms of resonant group of electrons. We 
reproduce this work in Chapter g. The new results of 
Chapter ^ are 

• resonant group of electrons which give the main 
contribution to the current, 

• "collisionless" current - analog of coUisionless Lan- 
dau damping of plasma waves. 

• difference between transverse magnetoresistance 
and magnetoconductance. 

The use of the Boltzmann kinetic equation in Chap- 
ter H has to be justified. In certain cases one has 
to use the more general quantum kinetic—squation as 
it was done by_Levinson and YasevichuteEf r-3uris and 
ShchamkhalovaEj and Laikhtman and Milleiu3. We will 
follow the last work in order to show that 
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FIG. 1.4. Vertical current through superlattice in presence 
of longitudinal magnetic field, (a) Current goes around thick 
areas of the barriers and becomes sensitive to the magnetic 
field.(b) 



• electron states in the narrow miniband in the pres- 
ence of the parabolic potential created by a mag- 
netic field are of two types Landau-like and Stark- 
hke, 

• the quantum kinetic equation supports qualita- 
tively results of the semiclassical approach of Chap- 
ter g 

• heating of the electron gas near the current peak 
can be important for wide interval of magnetic 
fields. 

The solution to the kinetic equation in the miniband 
is possible only if we are in the regime of the miniband 
transport. Four conditions of the miniband transport are 
given in the next section. It may happen, that some of 
these conditions are violated and we have to solve the 
quantum kinetic equation. In this case we have to justify 
also the effective Hamiltonian method and that is done 
in Sec. 



1.4 



The quality of the superlattices is determined basi- 
cally by the quality of the interfaces. Surface rough- 
ness destroys the coherence of electron tunneling, and 
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FIG. 1.5. Structure of energy levels in superlattice. (a) - 
Zero coupling and zero electric field; (b) - Coupling of quan- 
tum wells lifts degeneracy; (c) - strong electric field breaks 
miniband to Stark ladder. 



measurement of roughness parameters is very impor- 
tant technologically. It was suggested recently, that sur- 
face roughness can be responsible for longitudinal mag- 
netoresistance of superlattice. E3 If the superlattice lay- 
ers are isotropic, there should not be classical longitudi- 
nal magnetoresistance, (see the geometry of the experi- 
ment in Fig. 1.4(a)). Magnetoresistance was expla ined by 
the roughness induced in-plane currents, see Fig. 1.4(b), 
which are sensitive to the magnetic field perpendicular to 
the layers. Thepxy of this effect was developed by Miller 
and LaikhtmarO. We reproduce this work in Chapter ^. 
The main result is that longitudinal magnetoresistance 
of superlattices is very sensitive to the correlation length 
of the roughness and allows calculation of this length. 



1.3. FOUR CONDITIONS OF MINIBAND 
TRANSPORT 

In this section we introduce the main parameters char- 
acterizing electron states and electron transport in a su- 
perlattice. An example of a superlattice having = 4 
periods is shown in Fig. 1.1. Bold horizontal lines mark 
the positions of the electron levels in the quantum wells 
as if these wells are isolated. Then, the lowest states in 
all four wells have the same energy E = 0, and we mark 
it by single bold line in Fig. 1.5(a). 

Coupling of quantum wells lifts this degeneracy. The 
tunneling energy. A, depends strongly on the height of 
the potential barrier separating adjacent wells and on 
the effective mass of an electron under the barrier. Cou- 
pling of quantum wells spreads energies of electron states 
in the interval of width 4A, see Fig. p_.5K b). In the limit 
of the infinitely long superlattice N oo, this inter- 
val is covered densely by the levels. The energy interval 
[— 2A, 2A] is called a miniband. Therefore, by definition, 
a miniband exists in an infinitely long idealized super- 



{Ef,T)tp 



sequential 
tunneling 



miniband 
transport 



hopping 

conductivity 



Ar„ 



FIG. 1.6. Mechanisms, which govern vertical transport, are 
shown in the diagram. 



lattice, without impurities, interface roughness, phonons 
and external fields. 

An electric field, F, applied in the vertical direction, 
see Fig. L^, causes a relative potential shift of quantum 
wells. The drop of the electron potential energy per pe- 
riod is eFl — eV/N, where e is the electron charge, I 
is the superlattice period, N is the number of periods, 
and V is the applied voltage. If the coupling energy is 
zero, then this potential drop destroys the degeneracy 
of electron states and makes them equally spaced with 
the interval eFl, see Fig. |l.5| (c). This structure of the 
electron levels is called a Stark ladder. 

However, if both energies eFl and A are not zero, one 
can observe the transition from the leve l structure of 
Fig. ^(b) to the level structure of Fig. 1.5(c). This tran- 
sition is called the Stark effect. It has been intepapd* 
studied during the whole superlattice history.E3o'c2l 
Stark states have a finite size in the vertical direction 
Ay^A/{eF). 

A Stark ladder appears when the electric field satisfies 
the condition 



zFl > 



A 

N 



(1.1) 



i.e. when the mean spacing of the Stark ladder exceeds 
the mean spacing of the miniband. We see, therefore, 
that in the infinitely long superlattice the infinitesimal 
field destroys the miniband. Nevertheless, the picture 
of miniband transport is justified under the conditions 
which will be discussed in the rest of this section. 

The scattering of the electrons by charged impurities 
and by phonons can be characterized by the elastic Tp and 
the inelastic relaxation times. The scattering destroys 
plane waves of the free electrons and also Bloch waves 
of the electrons in the superlattice. However, the plane 
waves approximate the electron wave function well if the 
mean kinetic energy of the electron gas is large enough 



Ef,t:$> h/Tp,h/T, 



(1.2) 
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This is the first condition of miniband transport. If it is 
violated one can talk only about hopping between spa- 
tially localized states, which leads to hopping conductiv- 
ity. 



The condition Eq. (1.2) compares the kinetic energy of 
the electrons with the energy uncertainty due to scatter- 
ing. However, it may happen that the miniband width is 
smaller than this energy uncertainty. In this case the sub- 
sequent tunneling events lose their coherence and there 
are no more plane waves propagating in the vertical di- 
rection. Therefore, we obtain the second condition of 
miniband transport 



A > h/Tp,h/T^ 



(1.3) 



If this condition is violated then electron propagation in 
the vertical direction is called sequ ential tun neling. 

The two inequalit ies Eqs. (1.2) and (1^) define the 
phase diagram Fig. 



1.6 



In this diagram we show the 
three main regimes of super lattice transport. In the 
regime of miniband transport one can safely make use 
of the classical kinetic equation. In the regime of sequen- 
tial tunneling there are no plane waves propagating in 
the vertical direction and therefore one should compute 
transport coefficients from the quantum kinetic equation. 
Various percolation type models can be used for investi- 
gation of transport in the hopping regime. 

The strength of the electric field is also limited in 
the regime of miniband transport. It is clear that the 
Stark effect can be observed if eFl is much larger than 
h/Tp, Ti/t^, and-.that is very well known from the optical 
measurementfH. In the opposite case we can consider 
the electric field as a force driving electrons and that is 
precisely what we want. Therefore, we obtain the third 
condition of miniband transport 



^ h h 



(1.4) 



Violation of this condition leads to the so-called field- 
induced localization of the electrons in the Stark states 
and to the transport of the sequential hopping type. 
This raeans the hopping between the subsequent Stark 
states.El This regime cannot be observed experimentally 
because of the potential profile instability. 

The field induced localization is an interesting physical 
phenomenon. It has, of course, something to do with a 
Stark transition, but we would like to discuss it also in 
different terms. The main purpose of the following dis- 
cussion is to obtain the condition of miniband transport 
in the presence of a transverse magnetic field. 

Let us assume that all conditions of the miniband 



transport Eqs. (1.2), (1.3), and (O) are fulfilled. In this 



case one can introduce the semiclassical picture of the 
miniband inclined by an electric field. The energy of the 
electron having zero in-plane momenta Px = Pz = lies 
in the interval 




FIG. 1.7. Electrons reflects many times from the miniband 
boundaries before being scattered. The electron orbits are 
shown by thick lines (b). The Elastic scattering is accompa- 
nied by the large momentum transfer Ap (a), and therefore 
the inelastic process may be important. The inelastic transi- 
tion is shown by the dashed lines, see (a) and (b). 



eFy -2A< E < eFy + 2A 



(1.5) 
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This equation is semiclassically correct, because it con- 
tains both the coordinate y of the electron and the kinetic 
energy of the vertical motion of the elec tron . 

We show an inclined miniband in Fig. 1.7 on the back- 
ground of the supcrlattice potential. The relative shift 
of the subsequent quantum wells is eFl. In semiclassical 
language, the electron sitting at the bottom of the mini- 
band is accelerated by the electric field. It crosses the 
minibnad having passed the distance Ay « AA/{eF) and 
hits the upper boundary of the miniband. After reflection 
it goes back to the bottom of the miniband. 

The period of such a trip can be estimated by dividing 
Ay by the mean electron velocity in the miniband Al/h, 
which gives ~ T%/{eFl). The periodic motion of the elec- 
tron between the bottom and the top of the miniband is 
called a Bloch oscillator. It has tsftti observed in opti- 
cal four- wave mixing experiments J_3o The semicla ssica l 
language gives another physical meaning to Eq. (1.4), 



namely that the electron does not swing a long time be- 
tween the miniband boundaries because of the scattering. 

The swinging of the electron means that there is field 
induced localization, because the electric field is applied 
but cannot drive the electric current. This picture is very 
clear and there is no reason not to use the semiclassical 
theory of the miniband transport at least near the tran- 
sition point 



eFlTr, 



(1.6) 



Far from the transition point, where eFWp 3> h the main 
contribution to the current is given by the hopping be- 
tween subsequent Stark levels, and the contribution of 
the next neighbor hopping contains the small parameter 
h/leFl r). T his is not true at the instability point given 
by Eq. (1.6), where hopping between all Stark levels con- 
tributes to the current. A number of works, which ignore 
summation over all Stark levels near the instability point, 
have been published and shpw a wrong condition for in- 
stability and current peakc2l. 

It is known that elastic relaxation processes are usu- 
ally faster than inelastic ones. For this reason we put 



Tp in Eq. (1.6) and not r^. However, inelastic scatter- 



ing processes can be importa nt n ear the field induced 
localization transition. In Fig. 1.7(a) we showed the in- 
plane dispersion of the electron energy for two subsequent 
Stark states. The electron may need to transfer a large 
in-plane momentum to the impurity in order to jump 
between these states elastically. At the same time the 
inelastic process may be very effective. 

The acceleration of the electron by the electric field is 
equivalent to the presence of a linear potential, and to 
the motion of the electron in the tilted miniband, as it 
is shown in Fig. 1.7. The transver se m agnetic field cre- 
ates a parabolic potential, see Fig. l.S. In this case, the 
size of the electron orbit depends very much on the in- 
plane electron characteristic velocity (let the magnetic 
field be z-directed see Fig. 1.1(a) ). In order to see how 
it happens, we introduce the so-called effective Hamil- 




(a) 




(b) 



FIG. 1.8. The semiclassical picture of the superlattice mini- 
band bended by the transverse magnetic field, (a) - Bloch 
oscillator like orbits, (b) - Usual cyclotron orbits. 
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tonian; see the next section for the justificatfon of this 
approach. 



2m 2m 

Px ^ Px + eBy , 



2A 



Pyl 



eFy 



(1.7) 



where we used the MKS units and the tight-binding 
approximation for the dispersion law of the superlat- 
tice min iband. The coordinates here correspond to the 
Fig. [T^. 

Two integrals of motion Px and guiding center of the 
orbit yo allow us to introduce another conserving quan- 
tity, the in-plane characteristic velocity Vx = {Px + 
eByoj/m. We can obtain the one-dimensional Hamil- 
tonian for the particle with defined values of Vx and y^: 



n' = 2A 



1 



Pyl 



+ [eVxB - eF){y ~ yo) + 



[eB{y~yo)y 
2m 



(1.8) 



This equation shows that the kinetic energy of the verti- 
cal motion is limited by two inequalities 



E <2A+{eVxB~eF){y-yo)-\ 
E>-2A+{eVxB-eF){y-ya) 



°B{y~yo)]'^ 
2m 

[eBjy-yo)]' 

2m 




, (1.10) 



This 



which are quite analogous to the case of Ec 
energy lies between two parabolas, see Fig. 

There is a basic difference between two parallel lines 
and two "parallel parabolas" . Whereas the energy dis- 
tance between them is constant 4A, the size of the orbit. 
Ay, in the real space depends very much on Vx. The 
expression for the size of the orbit at the origin of the 
parabolas Fig. ^]^(b) is differ ent from the expression for 
the orbit in the branches Fig. 1.8(a). Therefore, we have 
two cases: 

\eF-ixB\ l^/^-^^l>\/S 
A2/=< • (1.11) 



\/2A^ 



\F/B-Vx\ < 



The first case corresponds to a slightly modified Bloch 
oscillator and the second case corresponds to the usual 
anisotropic cyclotron orbits. 

The periods of these orbits can be obtained by the divi- 
sion of Ay by the mean electron velocity in the miniband 
~ Al/h. The magnetic and electric fields induce local- 
ization of the electrons when these periods for most of 
electrons are of the order of the relaxation time. This 
happens either when 



^BTp 1 , hflB > \eFl - evpBll , 



(1.12) 



or when 

\eFl - evpBllTp n , TiVLb < \eFl - evFBl\ . (1.13) 

Here we replaced Vx by Fermi velocity vp- In the non- 
degenerat e cas e one has to substitute here the thermal 
velocity \fmT . The cyclotron frequency introduced here 
is defined by 



eB _ eBl A 
^mm± h V 2m 



(1.14) 



Note that the c ondit ion of the magnetic field induced 
localization Eq. (1.12) i s vah d for bulk semiconductors, 
while the condition Eq. ( 1.13|) i s spe cific f or sup erlattices. 

The analogy between Eqs. (1.12) and ( 1.13| ) from one 
side and Eq. (1.6) from other side leads to the fourth 
condition of the miniband transport: 

maxhriB, \eFl — eBvpl] i^nax — , — . (1-15) 



In order to provid e th e quantum mechanical arguments, 
as we did for Eq. ( |l.4|) , we have to quantize the effective 
Hamiltonian and compare the level spacings with relax- 
ation rates. The structure of the eigenenergies of the 
effective Hamiltonian is discussed in S ec. B. 13.1.1 and 
(1-9) one can check that conditions Eq. ( |1.15| ) are correct. 



1.4. EFFECTIVE HAMILTONIAN 

In the previous section we derived conditions of the 

(O), and (Osl). 



1.3) 



miniband transport Eqs. (|1^ 
The last one was obtained by making use of the effective 
Hamiltonian method. In this section we derive conditions 
justifying this method. 

The dynamics of a Bloch electron in the presence of 
a magnetic field can be described approximately-.by an 
effective Hamiltonian, introduced first by Peierlsa. The 
effective Hamiltonian is obtained from the spectrum of 
the Bloch electron i?(p) by the replacement of the mo- 
mentum p by p — eA, where|_A. is the vector potential of 
the magnetic field. Luttingeicijustified this approach for 
a weak magnetic field. BlountE3 calculated corrections to 
the effective Hamiltonian for a seperal three dimensional 
lattice. Berezhkovskii and SuriscZl made the accurate cal- 
culation of the electron spectrum in a superlattice and 
the external magnetic field applied perpendicular to the 
growth direction. They assumed that the effective Hamil- 
tonian can be applied if the magnetic field does not affect 
much the electron wave functions in separate wells of the 
superlattice. This results in the condition 



nVLcEp < El 



(1.16) 



where J7c = e.B/m, Ep is Fermi energy, and Eg is the 



energy gap between minibands, see Fig. 1.1(b). The con- 



dition of the applicability of the effective Hamiltonian 
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Eq. (1.16) can be considered as generally accepted, for 
the case of metals, see Rpf. 

Lifshitz and Pitaevskyca have assumed that the effec- 
tive Hamiltonian method can be used if the space be- 
tween energy levels in the presence of magnetic field is 
much smaller than the minimal energy scale of the band 
structure. We find a different condition. 



In this section we show that the condition Eq. ( 1.16 ) is 
not enough. We carry out the derivation of the effective 
Hamiltonian and find another necessary condition which 
can be more severe. We obtain our condition by the 
calculation of corrections to the effective Hamiltonian. 

In the absence of magnetic field, the eigenfunctions of 
the Hamiltonian Ho = {h'^\/^/2m) + U{y), where U{y) is 
a periodic potential, are 



1 



(1.17) 



where s is the miniband number and the Bloch amplitude 
Us is periodic in y. So, with the help of the unitary matrix 



k 



V 



V'.k(r') 



(1.18) 



the Hamiltonian can be transformed to the form where 
its wave functions are plane waves and the spectrum is 
the electron spectrum in the periodic potential, Es(k), 



Sl[^ + U{y)]Ss = EsHV) 



(1.19) 



In the presence of a constant vector potential this equa- 
tion becomes 



{hk - eA) 
2m 



+ U{y)]Ss^A^Es (-iV- Ja) 



(1.20) 



The matrix Ss,a is obtained from Eqs. ( |1.17| , |1 Jr^ ) by the 
replacement of Us{ky,y) by Us(ky — eAy/h, ?/)H3 

If t he v ector potential is not a constant then equation 
Eq. ( 1.20| ) is not exact. First of all, the matrix Ss.a 
introduced above, in general, is not unitary and has to 
be corrected. However, in a one dimensional superlattice 
and a uniform magnetic field it is possible to choose a 
gauge where Ay = 0, and the transformation matrix Sg 
is independent of the vector potential. Even with such a 
transformation matrix there are corrections to the right- 
hand side of Eq. ( 1.2C ) because Ss does not commute 
with the coordinate dependent vector potential. 

It's well known from the theory of Bloch wave functions 
that 



{s'k'\y\sky) 



(i^^+{s'\Y\s)yiky-k'y), (1.21) 



where (s|l^|s') = {us\i{d/dky)\us') is a function of ky. 
Here states {us\ are chosen in such a way that (s|y|s) = 



0. Thus, the diagonal matrix element of the Hamiltonian 
is 



2m 



(1.22) 



where B is the component of the magnetic field per- 
pendicular to the superlattice axis. The first term on 
the right-hand side of this equation is called the effec- 
tive Hamiltonian. In order to estimate the second term 
we make use of an equation for the operator Y. This 
equation can be obtained by calculating the commutator 
[[Ho, y],y] directly and with the help of Eq.(1.21), 



J2iEs + Es, - 2Es.)YssnYs»s' + i{Es - Es.)—Y, 



d_ 

dk' 



+2iYss>-{Es ~ E,) + (- - ^)<5..- = (1.23) 



If there are no band crossings, that is Ei < E2 < E3 < 
. . . for any fc, the diagonal part of this equation gives the 
following inequality 



11 



< 



l/m — l/mj_ 

T E2-E1 



(1.24) 



where l/m±{k) = d^Ei/{dk^). In the case of the smooth 
potential U{y) the right-hand side of this inequality can 
be considered as estimation for the mat rix e lement (i^^)ii 
in the left-hand side, see Appendix Sec 1.5. The variance 
of this matrix element in the Brillouin zone is, there- 
fore, of the order of {Ti^ /m){A/ Eg), i. e. it's propor- 
tional to the miniband width divided by the square of 
the energy gap between the first and second minibands. 
Eg =mmiE2- El). 

The last term in Eq. (1.22) in the case s = 1 can be 



neglected if it is much smaller than the miniband width, 
which gives the condition for the magnetic field 



hQc < Eg 



(1.25) 



Note that this condition is stronger than Eq. ( 1.16| ) unless 
the Fermi energy is comparable or larger than the Eg. 

The offdiagonal with respect to the miniband matrix 
element of the effective Hamiltonian can be neglected un- 
der the condition that is similar to Eq. (LIE). 



1.5. 



APPENDIX: MATRIX ELEMENTS OF Y IN 
THE TIGHT AND WEAK BINDING 
APPROXIMATIONS. 



In the tig ht bin ding approximation, A <^ Eg, both 



sides of Eq. ( 1.24 ) can be represented as a power series 
of the parameter 1/0 ~ A/ Eg. For example for a su- 
perlattice with rectangular wells, narrow barriers, and 
period tt 



1.5. 



APPENDIX: MATRIX ELEMENTS OF Y IN THE TIGHT AND WEAK BINDING APPROXIMATIONS. 14 



12 



TT 2 
6 TT 



2tt 



— + -]sm'{kn/2) +0{-^] (1.26) 



1 



For the same case the right-hand side of Eq. ( 1.24 ) is 



TT 

3 



1 1 

3 + e 

T ^ 97r / 



4 

97r 



20 \ 

sin2(fc7r/2)] +0 



1 

-Q2 



(1.27) 



So, in this particular case the calculated coefficients are 
different by about 3%. 

In the weak binding approximation, when band width 
is much larger than the gap, let 



1 



Un = - / cos(2ny)f/(y) , 



(1.28) 



where the superlattice potential U{y) is an even function 
with period tt, and Uq = 0. We assume that this potential 
is smooth; in other words we assume that the amplitude 
of the first harmonic t/i is large r tha n all other harmonics. 
In the gap region |fc| ~ 1 Eq. (1.24) becomes equality in 
the zero order perturbation theory: 



111 



1 



2 (1- |fc|)2+[/2 



(1.29) 



where we put h = m — 1. For small and intermediate k 
the second order perturbation theory gives 



oo 



1 



{k)~Es{k)dk^ \Ei{k) - E,{k) ^ 

(1.30) 



For the same case the right-hand side of Eq. (1.24) is 



dk^ \Ei{k)~Es{k) 



(1.31) 



Under the assumption that Ui > \Un\ for any n, the 
difi'eren ce be tween the left hand side and right hand side 
of Eq. drU) is less than 15%. 



Chapter 2 

Transverse magnetoresistance. I 



The motion of the Bloch electron in the crossed elec- 
tric, F, and magnetic, B, fields is studied for the case of 
the anisotropic band structure that usually exists in su- 
perlattices. The electron trajectories are calculated from 
the so-called effective Hamiltonian with an electric field 
applied in the growth direction of a superlattice and a 
magnetic field parallel to layers. We solve the kinetic 
equation and calculate the electric current for the case 
when the miniband width is much smaller than the char- 
acteristic kinetic energy of electrons. If the magnetic field 
is strong enough, only electrons that have a velocity com- 
ponent perpendicular to both electric and magnetic fields 
equal to the drift velocity, va = F/ B, contribute into the 
current. The current reaches its maximum value when va 
is close to the Fermi velocity, vp, or to the thermal veloc- 
ity, vt, whatever is larger. However, as the magnetic field 
goes to zero the current peak position goes to its limit 
Fth like B^ . The magnetoresistance at low magnetic field 
changes its sign at the field Fth/VS- The quantitative 
agreement with experiment is obtained for these results 
without fitting parameters. We also find that electric 
current is independent on the relaxation time in some in- 
terval of the applied field. This can be called the effect 
of coUisionless conductivity. 



2.1. INTRODUCTION 

The nonlinear transport phenomena in superlattices 
have attracted much interest in recent years. These phe- 
nomena come from the narrow width of electronic mini- 
bands in superlatticesa. Experiments show that with the 
increase of the applied bias the elecjtdc current first grows 
linearly, then reaches a maximun£3. At this point the 
uniform potential distribution across the superlattice be- 
comes unstable, and the superlattice breaks into low- and 
high-field domaina3. Further increase of the applied volt- 
age leads to the expansion of the high-field domain, which 
appears as a series of peaks on the I-V characteristic. 

Measurements of the vertical currenLj vi th e, p j res j en j se pf 
a magnetic field parallel to the layersH'BI^OSE3CZl'E3 
reveal a new interesting effect. It was found that all peaks 
on the I-V characteristic of a superlattice are shifted by 
magnetic field toward the higher bias. A linear theory of 



the magneto-conductivity of superlattices (i.<J»-| a theopy; 
for the Ohm's law region) developed by ShikEa, Andot2l 
and a theory of electron transport in rtiie Harper band 
developed by Suris and ShchamkhalovaEl cannot explain 
the behavior of current peaks. A qualitative explanation 
of this phenomenon was suggested by MovagjhaiEa. Nu- 
mericj-ealculations performed bpSibille et al.a. Palmier 
et al.E3, and Hutchinson et alH also lead to the same 
qualitative j^havior and some of them can be fitted to 
experimentEj. An analytical solution of the kinetic equa- 
tion and results for magneto-conductivity of superlattice 
have been publishedtHI. The final results were expanded 
in powers of magnetic field and they are reproduced by 
Eq. ( 2.25| ) here. In the present paper we calculate ana- 
lytically the I-V characteristic of a superlattice in strong 
electric and magnetic fields. The theory agrees quantita- 
tively with the behavior of the current peaks measured 
in different experiments. ■— ■ 

The qualitative arguments suggested by MovagharE3 
are based on the quantum mechanical explanation of 
negative differential conductivity (NDC) in superlattices. 
NDC in an electric field applied parallel to the superlat- 
tice axis {y axis) comes about because this field localizes 
electrons at Stark levels. A magnetic field applied paral- 
lel to the layers (in the z direction) causes the localization 
of electrons in the {x, y) plane. In a weak electric field the 
spatial size of electron states is limited by the magnetic 
field and it is smaller than the length of possible Stark 
states, A/ (eF). Therefore the electric field does not affect 
the electron states, and usual positive magnetoresistance 
is observed. In other words the I-V characteristic in a 
magnetic field goes below that without a magnetic field. 
With the increase of the electric field the size of Stark 
states eventually becomes smaller than the length of the 
possible magnetic states. At this electric field the effect 
of the magnetic field can be neglected and the I-V char- 
acteristic passes to the falling branch of the I-V curve 
without magnetic field. Because of the positive magne- 
toresistance, the peak of the I-V characteristic is shifted 
to a higher eleistric field. 

Sibille et alO calculated the I-V characteristic mak- 
ing use of the hydcpdynamical model similar to that used 
by Esaki and TsuB. The classical equation of motion for 
electrons in a superlattice in electric and magnetic fields 
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has been solved numerically. The electron drift veloc- 
ity was calculated by averaging of the electron velocity 
along the trajectory with some relaxation time. Although 
the hydrodynamic model gives a qualitatively correct I- 
V characteristic it docs not describe the experimental 
results, quantitatively. In a more recent work of Palmier 
et al£3 the I-V characteristic was obtained from the nu- 
meric solution of the Boltzmann kinetic equation in good 
agreement with the experiment. 

In this work we make an analytic calculation of the 
current in a superlattice in the presence of electric and 
magnetic fields with the help of the classical Boltzmann 
equation. The equations of electron motion in the crossed 
electric and n iagn etic fields in the superlattice are ob- 
tained in Sec. ^.2| from the effective Hamiltonian. These 
equations are siiialar to the equations of motion of a non- 
linear pendulurrO, but in our case two degrees of freedom 
are mixed and this complicates the solution. Electron 
trajectories are expressed in terms of Jacobian elliptic 
functions. In the Vx-Vy plane trajectories are closed or- 
bits confined in a strip \vy\ < 2M/h due to reflection 
from miniband boundaries. Here A is the tunneling ma- 
trix element between electron wave functions in adjacent 
quantum wells and / is the superlattice period. The cen- 
ters of the trajectories are located at the Vx axis. The 
size of each trajectory in the Vx direction is also limited 
from above, V&v{vx) < 4-y/2A/m, where m is the in-plane 
effective mass. The frequency of electron motion along 
an orbit has a minimum, when the center of the orbit is 
located at Vx = va = F / B. The frequency increases with 
the distance of the center from this point. 

In Sec. 2.3 we solve the Boltzmann equation in a su- 
perlattice miniband to calculate the electron distribution 
function, /. We assume that the characteristic electron 
energy (temperature or the Fermi energy) is much larger 
than the miniband width and the main scattering mecha- 
nism is elastic scattering by impurities or surface rough- 
ness. We do not make any assumption concerning the 
relative importance of external fields and collisions for 
the shape of the distribution function. This means that, 
in general, electrons may travel rather a long distance 
in physical space as well as in momentum space before 
being scattered. To solve the Boltzmann equation in this 
case we introduce as coordinates two integrals of mo- 
tion characterizing a trajectory in the momentum space 
and a coordinate characterizing electron position along 
the trajectory. The electron distribution is asymmetric 
along the trajectory with respect to the direction of the 
electric field. The asymmetric part is proportional to the 
relaxation time r when it is smaller than the period of 
motion along the trajectory and is proportional to l/r in 
the opposite case. 

The con- 



The electric current is calculated in Sec. 2.4 



tribution of an electron to the current is proportional to 
its velocity. The velocity in the growth direction, Wy, 
is an oscillating function of the electron momentum in 
the same direction, py. This momentum is changing in 
time under the electric field and the Lorentz force. For 



that reason Vy^ oscillates and the average contribution 
to the current is zero. However, the Lorentz force in y 
direction is proportional to the velocity in the direction 
perpendicular to the electric and magnetic fields, Vx^ and 
for some value of this velocity the force due to the elec- 
tric field and the Lorentz force cancel each other. This 
means that there is a resonant group of electrons with 
a time independent Py that contributes to the current. 
Scattering complicates this picture because for a nonzero 
contribution to the current it is enough that an electron 
does not change Vy between two scattering events. The 
width of the resonant region in the momentum space is 
determined by competition between the electron motion 
in the external fields and scattering. For some values of 
the fields a resonance region appears at the exponential 
tail of the electron distribution function. In this case a 
nonresonant contribution to the current, which is propor- 
tional to the probability for an electron to be scattered 
during one period of the orbital motion has to be taken 
into account. 

Application of a magnetic field perpendicular to the 
electric current generates an electric field parallel to the 
layers (in the x-direction) . This field is introduced into 
the kinetic equation in Sec. p.5[ We show that the cor- 
rection to the current in the superlattice direction due to 
this field (Hall effect) can be neglected. 

In a certain range of the magnetic field we found that 
the electric current is independent of the scattering rate. 
This result suggests that the kinetic equation can be 
solved without considering the collision term, as is done 
for plasma, see Ref. |6^. That is the regime of the col- 
lisionless conductivity considered specifically in Sec. 2.6. 
In Sec. 2.7 we compare the theoretical predictions with 



available experimental data. The summary of the results 
is given in Sec. 



2.8 



2.2. ELECTRON DYNAMICS 

We consider a superlattice in the tight binding approx- 
imation when, without external fields, the electron spec- 
trum in the first miniband is 

n2 



Pz 



2m 



2A 



1 — cos 



Pyl_ 

h 



(2.1) 



where m is the electron effective mass, A is the overlap 
integral between adjacent wells, and I is the superlat- 
tice period. In an electric field in the y direction and a 
magnetic field in z direction it is convenient to choose 
a gauge where the vector potential depends only on y, 
A = (-By, 0,0). 

The external fields can be approximately taken in to 
account by theijeffective Hamiltonian method introduced 
first by Peierlsuj. The effective Hamiltonian is obtained 
from the spectrum of the Bloch electron Ep by the re- 
placement of the momentum p by p — eA, (see justifica- 
tion of this approach in Sec. 1.4). The Hamiltonian has 
the form 
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Py 




Px 




FIG. 2.1. The trajectories in p- and in u-space. The quantity py ranges from —TiTi/d to -wh/d, px ranges from mF/ B — 18pA 
to mF/ B -\- 18pA, Vy ranges from —2v± to 2v±, and Vx ranges from F/B — Qp\/m to F/B + QpA/m. Trajectories in the p space 
are closed for < fig and open in the opposite case. The values of the ratio Q./Q,b for the trajectories shown at the figure are 
0.2, 0.5, 0.95, 1.0, 1.05, 2, and 3. 



Zm Zm 



+ 2A 



1 — cos 



Py 



eFy , 



(2.2) 



and gives the following expressions for the velocity com- 
ponents, 



the case ft > ^Ib there are two solutions for each value 
of fi. We have for first solution px < mF/ B and 



2n ,^ , 
Py{t) — — am(itt) 



Px{t) 



mF 



(2.7a) 
(2.7b) 



Px 



2v I sin ( ^ 



P 



- , (2.3) 



where Px = Px + ^By^ and — Al/h. The equations of 
motion resulting from the Hamiltonian Eq. (^.2[) are 



Px = eVyB = flnPA sin 



Pyl 



eB 



Py = eF - evxB = eF px 

" m 



where 



eB 



(2.4a) 
(2.4b) 

(2.5) 



PA = v2mA, and m^ — I (2AP) is the electron ef- 
fective mass at the bottom of the miniband. Equations 
(2.4) have the following integral 



sm 



{mF/B - Px) 



Pyl 

2h) ^ 4p2 



(2.6) 



where > is a constant of the integration. Equa- 
tions (2.4) are pendulum equations and their solution is 
expressed in terms of elliptic functions, see Fig. 2.1. In 



where the modulus of the elliptic functions is fc = U,b /^■ 
These equations describe open trajectories in momentum 
space. Since the values of the electron quasimomentum 
Py are limited by the interval —irh/l < py < nh/ l onl y 
a fractional part am(ili)/7r should be kept in Eq. ( 2.7a ). 
The second solution, which correspon ds to Px > mF/B, 
can be obtained from Eqs. ( 2. 7a , 2. 7b ) by changing fl to 

-n. 

In the case < Ub electron trajectories are closed and 
there is only one trajectory for each value of fl. Electron 
momenta Py{t) and Px{t) are again expressed in terms 
of elliptic functions, but the modulus is different, k = 

n/^B, 



, , 2a 
Py\t) — — arcsm 



sn{ilBt) 



mF i7 
Px{t) = — 2pA—cn{nBt) 

tS \Ib 



(2.7c) 
(2.7d) 



In the case Q, = Q,b the solution is expressed in the ele- 
mentary functions 



Py{t) — zb^arcsin [tanh(rit)] 

TTiF 

P-^i^) —5- T 2pAsech(r2i) . 

X3 



(2.7e) 
(2.7f) 
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For further calculations we also need expressions for 
the y-component of the electron velocity, 



4v±sn{flt)cii{flt) , 
Vy{t) = { 4v±ksn{QBt)dn{QBt) , 
4uxtanh(rit)sech(r2i) , 



k ^ ^B/n < 1 , 

k = n/ns < 1 , 

Q^Qb ■ 



(2.8) 



The motion of the electron in real space can be obtained 
by the inte gration of electron velocities. The integration 
of Eq. (2_^) shows that the y-coordinate oscillates around 



some value yo, which can be considered as the third in- 
tegral of the motion after and fl. The total electron 
energy is conserved and can be represented as a function 
of the integrals of the motion 



E 



2m 2m \ B 



2PA 



eFyo . (2.9) 



2.3. SOLUTION OF THE BOLTZMANN 
EQUATION 

For the calculation of the current we make use of the 
Boltzmann equation for one electron distribution func- 
tion /(p) 



eBvy- h (eF - eBv^) tt" 

dp^ ' dp 



y 



5/ 
dt 



(2.10) 



coll 



The collision integral on the right hand side of this equa- 
tion includes elastic scattering mechanisms (impurity 
scattering, surface roughness scattering and, possibly, ab- 
sorption and immediate emission of an optical phonon) 
and inelastic ones (acoustic and optical phonons). Elas- 
tic relaxation time is typically about lO^^^sec. Con- 
cerning inelastic relaxation, we assume that both the 
Fermi energy and temperature are smaller than the op- 
tical phonon energy hflLo so that the energy relaxation 
due to optical phonon emission can be neglected. ESCJ The 
inelastic relaxation time due to acoustic phonon scatter- 
ing typically is about 10~^''sec or larger. This means 



that the main relaxation mechanism in Eq. (2.10) is elas- 
tic scattering. 

We will consider the elastic scattering of the electrons 
as point scattering, when the electron position and elec- 
tron kinetic energy, Ep, are conserved. This approxima- 
tion is justified if the electron wavelength is much larger 
than the length scale of the scatter. For electron motion 
in the growth direction the electron wavelength ranges 
from one superlattice period to the superlattice length. 
At the same time, scattering in the adjacent wells can be 
considered as uncorrelated. The electron in-plane wave- 
length depends on the electron kinetic energy, that is 
typically about or less than lOmeV, and leads to the 
wavelength ^ 5000A. The length scale of the in-plane 



elastic scattering is very short because in superlattices 
the spacer separating impurities from electrons in quan- 
tum wells is smaller than the superlattice period and the 
correlation length of the surface roughn ess typically is 
about 500A (see e.g. Ref. ^and Section 4.4). 

It is important to note that Ep is not an integral of 
motion and is not conserved along a trajectory. Elastic 
collisions lead to relaxation of the electron distribution 
to its average over the surface of constant kinetic energy. 



m = 



Jdpf{p)S{E-Ep) 
J dpS{E ~ Ep) 



(2.11) 



The locality of elastic scattering means also that the ma- 
trix element of the scattering is a constant, which allows 
us to describe the elastic collision operator with the help 
of the constant relaxation time, [/(-Ep) — /(p)]/t. 

The electron gas can be appreciably heated by a strong 
electric field rigr '--^ 1. In this case f{E) can be signif- 
icantly different from the equilibrium function. In this 
work we will elaborate two limiting cases when the dif- 
ference between f{E) and the Fermi function with some 
effective temperature T, /o(£'), is not important for the 
calculation of the current. In the first case, T ^ Ep, 
the electron gas is degenerate and a detailed structure of 
the electron distribution near the Fermi surface does not 
make any difference. In the second case, Ep <C T w Ts 
where T^ is the lattice temperature, electron heating can 
be neglected.tj It is worth noting that at strong magnetic 
fields, quantum corrections to the classical Boltzmann 
equation ( 2.10| ) can be neglected only at high enough 
electron temperaturefH 



nnB 

2^ 



< 1 



(2.12) 



It is convenient to transform Eq. (2.1C)-to other vari- 
ables connected with electron trajectoriesEJ. These vari- 
ables are the integrals of motion 17 and pz, and time tp 
necessary for an electron to move from some fixed point 
in momentum space to the point p. In the new variables 
the kinetic equation becomes 



df _ h{Ep) - /(p) 
dU T 



(2.13) 



The formal solution to this equation is /(p) — /o(-E-p) -I- 
/i(p) where 



/i(p) 



eFvy{t)f'a[Ep{t))e-^'--''^'^dt . (2.14) 



Here /q is the derivative of the distribution function. 
Terms depending on time in Ep are of the order of A. 
If these terms are neglecte d the n /q can be considered 
time independent and Eq. ( ^.14 ) becomes 



h{p)=-eFTQp 



dfo_ 
dEr. 



(2.15) 
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where 



Vy{t)e~'-*''-'^/^dt = (1 + Td/dtp)-\y . (2.16) 



is a periodic function of time. The Fourier series repre- 
sentation of Qp is given in Appendix |2. 



Corrections to Eq. ( 2.15| ) contain A"/o which for- 



mahy is of the order of (A/T)". The electric current, 
however, contains integrals of this function and after an 
integration by parts, /q""''^^ in the integrand is replaced 
by a quantity of the order of //,/ niax( £'j?, T )". There- 
fore, for the calculation of current, Eq. ( ^.15| ) is justified 
when 



A < niax(£'F,r) , 



(2.17) 



where half of the perio d Tp is give n by Eq. (2.63). One 
can verify that Eqs. ( 2.19[ ) and ( 2.20 ) have the same 
asymptotic behavior when 1/t fts ^ ^■ 

These two limits show that the perturbation of the 
distribution function due to an electric field has a res- 
onant behavior. The physi cal r easons of the resonance 
have been discussed in Sec. If the resonance region 
is wider than the region of closed trajectories the maxi- 
mum value of Qp is of the order of v±. If the relaxation 
time is very long, O^r » 1, the value of is reduced 
by the factor of (fisr)^ due to oscillations of Vy along 
the electron trajectory during the time r. 



2.4. CALCULATION OF THE ELECTRIC 
CURRENT 



which is assumed to be satisfied in this work. 

The factor Qp can be separated into two parts, Qp ~ 
Qp + Q^ where Qp is even and Q^ is odd with respect to 
Vy. Only the odd part con tributes to the current in the y 
direction. From Eqs. (2^) and (2.4) it is obvious that Vy 



is odd in time so that the parity with respect to Vy is the 
same as the parity with respect to time. For that reason 



Q^ is obtained from the second expression in Eq. (2. If) 
keeping only even derivatives of Vy. Therefore 



Vp ^ 



1 



(2.18) 



This function is easily calculated in two limiting cases. 
Far from the region of closed trajectories given by the 
condition fl 3> ^b, we have Vy — 2v±sm{2ilt). Then 
d'^sm{2nt)/dt'^ = ~{2n)^sm{2nt) and we obtain 



1 + {2nTY 



(2.19) 



Under the same condition, 2^l « ilB\mF/B — Px\/PA, 
which means that Q^ as a function of px has a resonance 
at Px = mF/ B . The width of the resonance is pa/ (f^ s t) . 
If this width is much larger than the width of the region 
of closed trajectories, pA, th at is ^ 1, the latter can 
be neglected and Eq. ( ^.19 ) can be used for all values of 

Px- 

The second limiting case is a strong magnetic field 
when r^B T 3> 1, i.e., the width of the resonance in 
Eq. ( |2.19 ) is smaller than the width of the regi on of 
closed trajectories. Then the unit in Eq. ( 2.18 ) can 
be neglected compared to the second derivative, so that 
— r^(9^Q^/9t^) — Vy. The integration is carried out with 
the help of the equation of motion, Eq. (2.4), and the re- 
sult is 



Ql = 



2v±_ 



h Tp 

Pyd 



(2.20) 



n<nB 



The part of the distribution function which depends 
only on Ep does not contribute to the electric current, 
and the current in the field direction is 



J 



'2e^FT 



dp dfo 
(2nh)^ dEp 



(2.21) 



Two factors in the integrand, rf/o /dEp and Q^, have max- 
ima and the value of the integral crucially depends on the 
relation between the widths of these maxima. The width 
of the Q^ maximum depends on the value of fisr. If 
this product is small then the width is p\/ (^Ibt). In 
the opposite case it is pa- The width of the distribution 
function derivative is of the order of pr = \/2mT. Three 
cases are possible, (fisr)^ <C [pk/ptY — when 
the region of closed trajectories is not important and the 
width of Q^ is much larger than the width of dfo /dEp, 
A/T {^IbtY 1 when the region of closed trajecto- 
ries also is not important but the width of Q^ is much 
smaller than the width of dfo /dEp, and ri^r ^ 1 when 
only the vicinity of closed trajectories contributes to the 
resonance. We consider these three cases separately. 



2.4.1. Weak magnetic field, (fisr)^ < A/T 



In this case Eq. ( p. 19 ) can be used for computation 
of Q^ on the right hand side of Eq. ( ^.21] ). After the 
integration with respect to Py one obtains 



3v 



Ae'^FlrA^ f dpxdp. 



{2nnf 



X 



dh 
dE„ 



1 



PA 



(2.22) 



where we introduced frequency of the Bloch oscillations 

eFl 



Qp = 



h 
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Under the condition in the title of this subsection the ab- 
solute value of the momentum in the integral is controlled 
by dfo/dEp and 



Jy 



where 



d(j) 



1 + {flFT — flBTpF cos 4>/pa) 



(2.23) 



Pa 



r^ = {v + 2f + 4nlT^P^ 
Pa 



(2.24) 



In the nondegenerate case the main contribution into the 
integral Eq. (2.22) is given by px ^ pt which is small 
compared to pa/ ^bt . If Px is neglected in the square 
brackets in Eq. ( 2.22| ) the result is independent of the 
magnetic field. So in this case the magnetic field intro- 
duces just a_small correction computed for the first time 
by EpshtcinEl, 



Jy 



fiT 1 + {P-FTf 



2 2 T 3(r2FT)2-l 



1 + ^rt' — 



2A [1 



(2.25) 



wher e n is the s heet concentration. Both of the results, 
Eqs. ( 2.23| ) and (2.25) describe an I-V characteristic that 
has a peak. At zero magnetic field the peak position is 
given by the condition VlpT = 1. With the increase of 
the magneti c field this peak shifts to a higher bias. 

Equation (2.25) shows that magnetoresistance changes 
its sign at nj.r'^ = 1/3, or F = hl/{eTV3). The mag- 
netoresistance is positive for a weaker electric field and 
negative for a stronger one. For the Fermi statistics the 
same effect takes place when the magnetic field goe s to 
zero, which can be proved by expansion of Eq. (2.23). 




mF/B 



FIG. 2.2. The division of the electron distribution into a 
nonresonant part and a resonant part (hatching area) in an 
intermediate and strong magnetic field. The width of the 
resonance, Ap, is pA or pa/{^bt), whichever is larger. 



^ J (mr/7r)3/2 

PA 



PAG 



^pj, - {mF/Bf 



, t:s>Ef 

T-^Ef 
(2.26b) 



In the case of Fermi statistics, Eq. ( |2.26bD is valid only if 
PF — mF/ B ^ Pt. This result is also the limiting form of 
Eq. (|l23|) when pp-mF / B :^ pa/ (S^rt). Therefore the 
reciprocal square root singularity obtained near the point 
mF/B = Pi? is smeared either due to the temperature or 
due to the width of the resonance, whichever is larger. 
When mF/B > pp the current drops off and the position 
of the current maximum is given by mF/B = pp, i.e., 
F = {pp/m)B. In the case of Boltzmann statistics, the 
position of the current maximum is mF/B = px/V^, 
i.e., F — {pT/m\/i)B, meaning that in the intermediate 
magnetic field limit the current peak shifts linearly with 
B. 



2.4.2. Intermediate field, A/T < [yiBrf < 1 



2.4.3. Strong magnetic field, Q.bt 3> 1 



In this case Eq. (2.22) also can be used. However, the 
value of Px is controlled now by the Lorentz factor. This 
factor picks up particles from the phase space, which have 
valu es o f momentum near px — pa^f/^b = mF/B, see 
Fig. |2.2|. Then the integration with respect to Px gives 



HF,B) , 



(2.26a) 



p,=mF/B 



The calculation in this case is similar to that in the 
previous section. The only difference is that it is neces- 
sary to use Eq. ( 2.20| ) for and the integration with 
respect to Py can be carried out in a different way. After 
the separation of the integration with respect to Px 



8R e^mA^l 



(2.27) 



where the dimensionless constant of the order of unity, 

„2 ^ -1 



R 



327r 



I {^BTf 



Qlvydpxdpy , (2.28) 
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is calculated in Appendix. 2.10| . 



It is instructive to note that Eq. ( 2.27 ) can be written 
down in the form similar to that in the bulk, 



F 



J 



(2.29) 



Here h is the effective electron concentration, i.e., the 
number of the electrons inside the resonance region 



^ HF B) 



(2.30) 



eBi 



JH 



J 



(0) 



(2.35) 



which is very well known for the bulk semiconductors 
with a constant relaxation time. 

In order to calculate axx and ayx it is necessary 
to introduce into the Boltzmann equation a small x- 
component of the electric field and calculate the pertur- 
bation of the distribution function, /2, caused by this 
field. After the linearization with respect to this field, 
the Boltzmann equation becomes 



dpx dtp 



A 

T 



(2.36) 



2.5. ESTIMATION OF THE ROLE OF THE 
HALL EFFECT. 

Typically, in the current-voltage characteristic mea- 
surement in superlattices, Hall contacts are not fabri- 
cated. Then, in a magnetic field, a Hall voltage comes 
about which can affect the value of the current. As a 
result of this effect in a weak electric field the quantity 
measured in the experiment is not the conductivity axx 
but the resistivity pxx- In this section we show that the 
effect of the Hall voltage on the current-voltage charac- 
teristic in a superlattice can be neglected. 

The Hall field Fx and the current in the growth direc- 
tion can be calculated from the equations 



Jy = Jf' + 

jx ^ jH + ' 







(2.31) 
(2.32) 



where jy^"^ is the value of the current for zero Hall field 
and Jh is the Hall current under the same condition. All 
the components of the conductivity may depend on the 
field Fy. Both jy^^ and jn nontrivially depend on Fy. 
The conductivity components dyx and a^x also can de- 
pend on Fy . This dependence results from a narrow mini- 
band characterizing the electron motion in this direction. 
There is no such small energy scale for the in-plane elec- 
tron motion, so that the effect of Fx on jy^\ ju, and the 
conductivity components can be neglected. 
In the calculation of the Hall current 



in = 2e 



dp 



(27r?i)^ 



(2.33) 



it is convenient to make use of the identity 



fi{p)dpy = -eBr— Vyfi{p)dp. 

-■nh/d OPx J-Trh/d 



7vh/d 
h/d 



(2.34) 



which is obtain ed by the integration of the Boltzmann 
equa tion ( 2.10|) wit h resp ect to Py. The substitution of 
Eq. (|2.34D into Eq. (|2.33D gives the resuh 



where / = /o + fi is the distribution function obtained 
in Sec. 2.3 for the case F^ = 0. 



The integration of Eq. (2.36) with respect to py gives 
the identity 



7Th / d 



lip)dPy 



— nh/d 



-eFxT- — 

OPx 



eBr— / vyh{p)dpy , (2.37) 

OPx J-nh/d 



that results in the relation between the conductivities 

eBr 



ne^T 



md 



(2.38) 



Therefore only one of them, say Uyx , has to be calculated. 
A formal solution to Eq. ( p. 36 ) is 



h 



-eFxT 



d 



rd/dtp dpx 1 + Td/dt^ 



-k{Ep) 



(2.39) 



Then 



d 



l + Td/dtp dpx l + Td/dtp 



UEp) 



where 



^2eV 



dp 



(2.40) 



(2.41) 



We will show now that tTyx can be expressed in terms of 

CTyy which is defined by the re lation j'^'^ 
the help of the definition Eq. (2.41) 



ayyFy. With 



'yy 



dfo\ 
-V,,—— > . 



1 + rd/dtp " dE I 
With the help of the commutation relation 



d_ 
'dL 



d 
dpx 



eB d 
m dpy 



(2.42) 



(2.43) 



it is easy to show that 
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l+T 



~ _d_ _ _d_ 

dpx dp: 
d d 



dpy 



_d_ 



(2.44) 



The substitution of Eq. ( |2.44| ) into Eq. ( |2.40| ) gives 
eBr 



d 



l + rd/dtpdpy \l + Td/dtp 



Now we make use of the commutation relation 



_d_ 



d 

dpy 



eB dvy d 
m dpy dpy 



(2.45) 



(2.46) 



to estimate the commutator of the operator factors in 



Eq. (2.45). If we take into account that d/dtp ~ fl then 



the ratio of the commutator ofd/dpy and {l + rd/dt^ 
to the product of these operators is of the order of 
[17bt/(1 + Jlr)] (pa/Px)- T he la st quantity is much 
smaller than unity due to Eq. ( |2.17| ) and 17b < $7, which 
means the commutat or of d/dpy and ( 1 + T d/dtr,)~^ can 
be neglected in Eq. ( 2.45 ). Then Eq. ( 2.45 ) becomes 



cBt 



1 



' dE 



(2.47) 



The expression above is proportional to the second 
derivative of the (Xyy with respect to the field, which can 



be proved as follows. First of all Eq. (2.42) gives 
52 




li^ ■ (2-48) 



Here we neglect the commutator of d/dpy and (1 + 
Td/dtp)~^ once again, and make use of d^Vy/dpy — 



Eq. (2.48) gives 



The comparison of Eq. ( 2.47 ) with 



eBr 
2m 



n 

erl 



yy 



dFy' 



(2.49) 



In the limit of a small Fy one can easy verify that — Cy 



CTx 



dn/Fy. 

Now we can estimate the value of the second term in 



Eq. (2.31). Making use of Sec. 2.4 we have the relation 
of the second term in Eq. (2.31) to the first one 



cTy^Fx _ eBr a 
jo r 



yx 



^ Ef,T ' 



A 



1 



A 



Qbt \Ef,T 
1 / A 



5/2 



5/2 



A 



Ef.T 



Ef,T 

< f^^jT^ < 1 



17|t2 \Ef.T 



1 < n^T^ 



(2.50) 



and one can see that this ratio is small in all considered 
limits. 



2.6. DISCUSSION 



One of the most interesting results of Sec. 2.4 is that 



in some region of magnetic field, A/T ^ {VIbT) <C 1 
the current does not depend on the relaxation time, see 
Eq. (2.26). That means that a finite resistivity of a su- 
perlattice exists even without scattering. To understand 
the physical reason of this result it is instructive to note 
that the considered problem is quite similar to the prop- 
agation of a longitudinal wave in a collisionless plasma. 
Indeed, if we consider the miniband width as a small 
parameter {vy = [2M/Ti\sva[pyl/Ti\ and for A = the 
current in the growth direction is zero) then in the linear 
approximation in A and neglecting collisions, Eq. ( 2.10| ) 
takes the form 



eB 



df_ _ 2eBM ^^^pyl_ dfojEp) 
dpy h h dpx 







(2.51) 



On the other hand, the Boltzmann equation for a col- 
lisionless plasma rwith an electric field in the wave 
Fq sin k{x — Vpt) isE3 



{vx-vp)^ + eFo sin k{x - Vpt) = q . (2.52) 

ox Opx 

The comparison of these equations shows that there is a 
one to one correspondence 



plasma 

X - 

k 



Vpt 



eFo 



superlattice 
Py/eB 
eBl/h 
F/B 

-2eBAl/h 



In a plasma there exists a resonant group of electrons 
moving with the velocity close to the wave velocity, Vp. 
This group of electrons strongly interacts with the wave, 
resulting in the collisionless Landau damping. A similar 
resonant group of electrons exists in a superlattice. Those 
are electrons moving in the x direction with the Hall 
drift velocity, F/B. For those electrons the Lorentz force 
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\B 



II 



Bth 



III 



A 



FIG. 2.3. I,-region of the strong magnetic field, II,-region 
of the colhsionless conductivity, and 111,-region of the weak 
magnetic field. The threshold value of the magnetic field, 
Bth ~ hm/[elT{pF,PT)], is only a few T for superlattices. 



cancels out with the force of the electric field and their 
velocity in the y direction does not oscillate in time. As 
a result their contribution to the current is nonzero. 

The correspondence between a monochromatic plasma 
wave and a superlattice in crossed electric and magnetic 
fields can be extended even beyond the linear Landau 
damping in plasma and the linear in A theory in a su- 
perlattice. Because of some specifics of collisions in a gas 
plasma, a close correspondence in this case exists between 
the superlattice and the solid state plasmaO. We show 
here only the correspondence between the applicability 
conditions for the linear theory. 

In plasma a wave with a finite amplitude traps the res- 
onant electrons and they oscillate in potential wells of the 
wave. The linear theory is applicable if the period of the 
oscillation is larger than a scattering time which can be 
written as eF^kr'^ /m <g; 1. The above correspondence 
between the plasma and superlattice quantities immedi- 
ately gives the condition for the superlattice, ^%t^ ^ 1. 
The physical meaning of that condition is quite similar 
to the plasma one. A finite width of the miniband in the 
superlattice leads to Bloch oscillations of resonant elec- 
trons. Then the theory which is linear in A is justified if 
the period of the oscillation is larger than the relaxation 
time. 

The situation when the collision term can be neglected 
in the kinetiCpSquation is typical for plasma and given by 
the conditionEJ kpTr/m ^ 1. Regarding the superlat- 
tice, this condition becomes the condition of the inter- 



mediate magnetic field ^bt{j>t /ph) ^ 1- In the case of 
the degenerate electron gas important for superlattices, 
the collisions ca n be n eglected under a more weak condi- 
tion. Equation. ( 2.22 ) takes the form of the Eq. (2.26) for 
Ep ^ T, when VIbt{pf /Ph) ^ 1, i-c the range of the 
magnetic field when the coUisionless conductivity takes 
place in superlattices is 



A 



Ef,T 



2 ^2 



< 1 



(2.53) 



The existence of this effect depends on the relation be- 
tween miniband width and characteristic electron energy 
which is schematically shown in the diagram Fig. 2.3. 
The plane B-A. in this diagram is divided into three parts 
corresponding to the different regimes. The coUisionless 
conductivity (intermediate magnetic field regime) can be 
obse rved only if A <C Ef,T, region II in the diagram 
Fig. |2.3| . The strong magnetic field limit, region I in 
Fig. ^.3^ is the same for small and large A; electric cur- 
rent exhibits 1/r dependence. The positive or negative 
magnetoresistance can be observed in the regime of the 
weak magnetic field, region III in Fig. 2.3, depending on 
the electric field. Due to the large period of the superlat- 
tice the boundary between regions II and III lies in the 
experimentally available range of the magnetic fields. 



2.7. COMPARISON WITH EXPERIMENT 

The influence of the electric and magnetic flelds on the 
current in the superlattice was measured in ±sKa, types of 
experiments. In the first kind of experimentsO'LD the I-V 
curve of the superlattice was measured in the presence of 
the constant magnetic field. The position of the current 
peak depends on the magnetic field and this dependence 
can be plotte^jm the B-F diagram. In the other kind of 
experimenta£3'E3 the dependence of the electric current 
on the magnetic field was measured in the presence of a 
constant electric field. A peak in the I-B curve has been 
detected, and its position can be plotted as a function of 
the electric field in the B-F diagram. 

We describe now the predictions of the present theory 
for such experiments. The peak on the I-V curve can be 
obtained from the equation dj/{dF) — 0. In the absence 
of a magnetic field the position of the peak is given by 
r^i^r = 1, or = Fth = h/ (elr). When a weak magnetic 
field is applied, B<^Bth, where 



B 



th 



hm 

epplT 
/— hm 



eprlr ' 



Ef->T 
Ef-^T 



(2.54) 



the peak position is shifted by small amount F — Fth <C 
Fth- In the lowest order of the mag netic field the shift 
obtained from Eqs. ( |2.23| ) and (|^) is 

^ eB^d ina^{EF,T) 
F- Ftb^i ^rn ■ (2-55) 



2m 



h/r 
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FIG. 2.4. The position of the maximum in I-V characteristics of the superlattice, empty circles, and on I-B characteristics, 
filled circles. The solid lines are theoretical predictions. 



With the increase of the magnetic field F — i^th becomes 
comparable with Fth- This happens when B^Bth and in 
this limit current is given by Eq. (2.26). In this interval 
of the magnetic field the peak position changes linearly 



F : 



(vf,vt/V2 



B = ^B 
Bth 



(2.56) 



The peak in the I-B characteristic is given by the 
equat ion dj/d B = 0. In the limit B goes to zero, 
Eqs. (2.23, 2.25) give the solution 



F 



(2.57) 



When B>Bt\i, the solution to equation dj/dB — co- 
incid es wi th the solution to the equation dj/dF — 0, 
Eq. ( ^.56 ), because the current is a function of F/B, 
see Eq. (2.26). That is true unless fisr > 1. In th e 
last case the current acquires the factor 1/B, Eq. ( p^ , 
and the peak in the I-B curve is shifted to weaker 
B. In the case of the nondegenerate electron gas and in 
the limit ilsr 3> 1 we have the peak position given by 
F = vtB. Therefore, the peak in the I ~ B curve lies 
between F/{^/2vt) and F/vt depending on SIbt. In the 
experiment the peak position in the I-B curve is mea- 
sured as a function of the applied bias, V = LF, where 
L is the superlattice length, and we have the following 
useful inequality for the bias corresponding to the peak 
position 



T 

'—LB <V < 
m 



—LB, B>Bu, 
m 



Ep^T . (2.58) 



In the work of Choi et aS the B - F diagram was re- 
ported. The superlattice investigated in this experiment 
has TV = 49 periods, the length of the period is I = 166A, 
the width of the first miniband is 4A = 0.47meV, and the 
concentration of the carriers corresponds to the Fermi en- 
ergy Ep = 6meV. At zero magnetic field the current peak 
appears near the bias Vth = NlFth — 590mV, and there- 
fore the momentum relaxation time is r = hN / (eVth) = 
0.54 X 10^^'^sec. The critical value of the magnetic field 
is Bth — 4.2T. The exp erim ental data from Fig. 7 of this 
work is plot ted in Fig. 2.4(a), empty circles. The solid 



lines on Fig. 2^ (a) are obtained from Eq. ( 2.23 ), since in 
this experiment Q bt <C 1 for the whole range of magnetic 
fields. The upper curve represents the maximum position 
on the I-B characteristic and the lower curve represents 
the maximum position on the I - F characteristic. One 
can see a good agreement between the lower curve and 
the experimental points without any fitting parameters. 

The both types of measurements (peak position on 
the I-B curve versus electric field and maximum of 
the I - F characteristic versus magnetic-field) have been 
reported in the work of Aristone et aEll. This experi- 
ment has been done on low doped samples with Ep much 
smaller than the temperature T = 23meV. The narrow- 
est miniband sample 4A = 4meV has a period / = 74A 
and a length Nl — l/im. The critical value of the applied 
bias is Vth = 500mV, and therefore the momentum relax- 
ation time is r = hN/{eVth) — 1-7 x 10~-^^sec. Thus, the 
critical value of the magnetic field is i?th = 2. IT. The ex- 



perimental data from this work is plotted on Fig. 2.4(b). 
The empty circles represent the peak position on the /- 
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V characteristic taken from Fig. 4 of Ref. |5^, and the 
filled circles represent the peak position on the I-B char- 
acteristic taken from Fig. 6 of Ref. On t he same 
graph, Fig. ^.4| (b), we draw the line Eq. ( 2.56 ) that is 
asymptotic to the peak positions for intermediate mag- 
netic field. The theoretical prediction for the weak mag- 
netic field regime is not shown on the graph, but one can 
see that the filled circles (the peak positions on the I-B 
characteristic measured experimentally) converge to the 
point B = 0, F = Fth/VS in agreement with Eq. ( 2.57 ). 
The good agreement between theory and experiment is 
obtained without any fitting parameters. 

In the experiment of Aristone the samples with the 
wider minibands were investigated too. The position of 
the peak on the I-B characteristics of these samples ex- 
hibits a weak dependence on the miniband width, that 
is in agreement with the present theory, since for these 
samples the limit Qbt ^ 1 is reached. One can check 



that the inequality Eq. ( 2.58 ) is fulfilled for this data. 



2.8. CONCLUSIONS 

In conclusion, we solved analytically the kinetic equa- 
tion for the superlattice in crossed electric and magnetic 
fields, under the assumption that the miniband width is 
small compared to the Fermi energy or the temperature. 
Our calculations of the position of the maximum on the 
I-V and I-B characteristics are in quantitative agree- 
ment with experiment. We also found that in a certain 
range of the magnetic field the current is independent of 
the electron relaxation time. 



2.9. APPENDIX A 

The Fourier series expansion of the velocity is 

2 oo 



Vy{t) = 4:V±_ 



^bTp 



E I TllnK'X . ( TTl„t 
t„sech I — — — I sm 



K 



p 

(2.59) 



The application of the identity 
1 



1 + Td/dt 



, , , sm(Lut) — LOT cosiut) „„x 
sm(.t) = (,,)2 (2-60) 



to each term in the series Eq. (2.59) results in the follow- 
ing series representation for the factor Qp: 



Av I 



^bTp 



2 oo 



Z„sech 



TTlnK'jk) 

K{k) 



sin (TT^ntp/Tp) TxlnT COS {Trl„tp/Tp) 



l + {TTlnTlTpY Tp i^i^^i^r/TpY 



n Q > flB 

n-1/2 n < Hb 



(2.61) 
(2.62) 
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max (Q, ^Ib) 
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^b 

n 



(2.63) 



Here K{k) is a Complete Elliptic Integral of the first kind, 
and K'(k) = K{y/1- P). The division of the factor Qp 
into parts odd and even with respect to Vy cor respo nds 
to the two terms in the square brackets in Eq. (2.61). If 



one keeps only the left term then the sum will result in 
Q^, and if one keeps only the right term then the sum 
will result in Qp. 
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The substitution of Eq. (2.20) into Eq. (2^28)) gives 



R = (167r ApaY 



Pyl ntp 2 



Vydpxdpy , (2.64) 



where 0{x) is the step function. The integral with re- 
spect to Py of the first term in the integrand is trivial 
but the next integration of this term with respect to px 
diverges. In the second term it is convenient to carry out 
the integration in variables 17 and t. The integral with 
respect to Q, also diverges. The sum of the integrals has 
to converge. So it is convenient to put the limits of the 
integration in the first term p^ = —pM and Px = Pm- 
In the second term the limits are = ili?/2 + f7M and 
Q. = Q.p/2 — D.M respectively, where VLm ~ ^bPm /'^Pa- 
The result is obtained in the limit pM oo. Then keep- 
ing in mind that the Jacobian of the transformation from 
variables , Py to il, t is 2{Q./VIb){piJi/1) and calculat- 
ing the integrals of the elliptic functions with the help of 
Eqs. (5.136.1) and (5.133.3) of Ref. ^we have 
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Pm 

2pA 



n%K{nB/^) 



n-2K{nB/n)\li-^ 



j2jj+riF/2 



n^dfi 



n%K{nBin) 



n -2K{nB/fl)\ 1- 
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(2.65) 



When n/rts goes to oo the integrand here goes to 1/2. 
This limit value can be integrated separately and then 



R=l 



1 + AkVk^ - 1 - 



27rfc2 
K{l/k) 



dk w 0.747 



(2.66) 



Chapter 3 

Transverse magnetoresistance. II 



In this chapter we derive the quantum kinetic equa- 
tion, which is a vahd approach for any relation between 
four energies A, Ti/Tp, eFl, H^Ib- Such a general case 
can be described quantitatively with the help of the den- 
sity matrix if an effective electron temperature is large 
enough. The qualitative results are similar to the results 
of Chap. H there exists a resonant group of electrons 
with Vx near F/B which gives the main contribution to 
the current. 

We showed in Chap. ^ that the current has a max- 
imum when F/B is close to the characteristic electron 
in-plane velocity. Therefore, the ratio F/B near the cur- 
rent peak allows one to measure such properties of the 
in-plane electron distribution function as electron tem- 
perature or Fermi energy. If the width of the miniband 
is smaller than the energy of opjtiftaJ phonons, the cool- 
ing of electron gas becomes slowE3lHa. The effective elec- 
tron temperature in the region of applied fields near the 
current peak can, be significantly larger than the crys- 
tal temperaturec3. Therefore, we need to have theoreti- 
cal predictions for this temperature; it can be compared 
with experimentally measured m[F / B^ near the current 
peak. 

Therefore, the purpose of this chapter is two-fold. We 
want to repeat part of the calculations of Chap. || starting 
from the quantum kinetic equation, and we want to see 
the effect of heating of the electron gas. 

Levinson and YasevichuteE^ considered for the first 
time the heating of an electron gas in an anisotropic semi- 
conductor or superlattice. They wrote down the quan- 
tum kinetic equation and considered simultaneously two 
effects: Stark localization and the heating of electron gas. 
Their results were obtained in the limit eFl ^ Ti/t, far 
away from the current peak; this region of fields cannot 
be observed experimentally in superlattices due to the in- 
stability of the potential distribution, see Fig. |1.2| and 
Chap. H The S|atne shortage is in the work of Suris and 
ShchamkhalovaEJ, who concentrated on the interplay be- 
tween the energy of optical phonons and miniband width. 

The calculations in|-t|his Chapter follow the work of 
Laikhtman and MilleiO. The presence of the magnetic 
field is taken into account in the appropriate places. Esti- 
mations for electron temperature and heating conditions 
near the current peak remained the same. 



We use the Keldysh technique in order to derive the ki- 
netic equation for the density matrix. We need the Fermi 
energy or the temperature to be much larger than the 
energy uncertainty F, that covers the upper part of the 
diagram in Fig. 1.6. In order to obtain analytical results 



we need to expand the kinetic equation with respect to 
A/Tg, and therefore the region of the miniband transport 
in Fig. 1.6 is covered only partially. Scattering is consid- 



ered in the Born approximation. We do not usj3 a trans- 
lationary and gauge invariant Green functionO, however 
we use a translationally and gauge invariant counterpart 
of the density matrix. 



3.1. QUANTUM KINETIC EQUATION. 

3.1.1. Quantization of tlie effective Hamiltonian. 

In this subsection we discuss the wave function of an 
electron moving in a one dimensional periodic superlat- 
tice potential U{y) in presence of the crossed electric and 
magnetic fields. For the purpose of this work we want to 
know how the external fields modify the wave functions 
of the lowest miniband produced by the superlattice po- 
tential. This problem can be solved by the meth od of 
effective Hamiltonian already discussed in Sec. 



1.4 



The method of effective Hamiltonian does not contain 
limitations on the shape of the potential U{y), and in this 
chapter we consider tight-binding case, sinsc it is most 
relevant to experiment. In order to quantize the effective 
Hamiltonian in this case, it is convenient to start calcu- 
lations in the Wannier representation, which is diagonal 
with respect to electron states in different potential wells 
of U{y). The explicit form of the effective Hamiltonian 
obtained in this way will be equivalent to that which was 
derived in Sec. 1.4. This can be proved by transformation 



to Bloch waves. 

Electron wave function in the Wannier representation 
is 5'-i/2g«Pi!z/?i-(-jPa:a:/?iy^(-y _ q,;-). ^]^jg jg ^j^g wave func- 
tion of the first level in the a-th well of U{y), I is the 
period of the superlattice potential U{y), (a;, z) are in- 
plane coordinates, S is normalization area. The electron 
Hamiltonian of a perfect superlattice in a uniform elec- 
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3.1. QUANTUM KINETIC EQUATION. 

E 




(a) 



(b) 



Landau levels 



FIG. 3.1. Two ways to quantize the effective Hamiltonian. 
See explanation in text. 



trie field F and z-directed magnetie field B, see Fig. 
pI ^ [P^+eBalf 



IS 



2 m 



eFla Saa' 



(3.1) 



Here e and m are the electron charge and mass re- 
spectively, and A is the tunneling matrix element be- 
tween first levels in nearest wells; we neglect tunneling 
on the distance more than one period. We neglect inter- 
miniband matrix elements of coordinate y, and presence 
of other bands; we also neglect field induced corrections 
to the tunneling. That is possible if the fields F and B 
are so small that eFl and eBvpl (wp-Fermi velocity) are 
much smaller than the distance between energy levels in 
a well. I— I 
The traditional way to quantize this Hamiltonian iscZI 
replacing the operator y or a l by —ihd/dpy. The field 
part of Hamiltonian Eq. (3.1) becomes a parabolic "ki- 
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netic" term proportional to —d^/dpy and the tunnel- 
ing term of Eq. ( |3.l| ) becomes a periodic "potential" , 
Asm{pyl/h). It is shown in Fig. |]^(a), that quantiza- 
tion of such Hamiltonian leads to formation of Landau 
levels at the parabolic bottom of each Brillouin zone. 
However, the tunneling between neighbor Brillouin zones 
results in the broadening of Landau levels into so-called 
Harper bands. These Landau levels type states are im- 
portant for transport if the Fermi energy or temperature 
of the electron gas are much smaller than the miniband 
width A. Conductivity of the lowest Harper band was 
computed by Suris and ShchamkhalovaEil by making use 
of the Kubo formula. 

In the opposite case of the large mean kinetic energy of 
electrons, the greatest contribution to transport is given 
by the states far above the periodic "potential" shown 
in Fig. 3T(a). These states are not similar to Landau 
levels at all. Berezhkovskii and Suris compute them from 
a weak binding model for "potential" Asm{pyl/h). We 
will use another approach in order to give new physical 
interpretation to these states. 

In the framework of the semi-classical approach, one 
can draw in the same picture miniband boundaries and 
the potential of an external field bending this miniband, 
see Fig. |3.l| (b). We immediately resolve in this picture 
two types of states: Landau-like states and Stark-like 
states. Near the bottom of the low parabola, the Landau- 
like states are formed (states of the harmonic oscillator) . 
They are slightly broadened by the presence of the top 
of the miniband, nevertheless they are shown by parallel 
horizontal lines with an energy spacing hflB- Far away 
from the center of the parabolas, the two parabolas are 
almost par allel near some particular energy, see dashed 
lines in Fig. |3.l| (b) . Electron wavefunctions are Stark- like 
near such an energy. Energy separation between Stark- 
like states is the slope of these parabolas multiplied by 
the superlattice period. 

A Stark-like state is localized around a certain well. 
The position or number of such a well is a good quantum 
number, and let's call it the guiding center of the elec- 
tron orbit v. After adding and extracting a few terms 
containing i', the Hamiltonian Eq. (|]^) can be rewritten 
in the form 



Ho 



2 2 

Zm 2 



(3.2) 



n'o - ( "^^^ + [eBvJ - eFl] {a - j,) ] 5^ 



2m 

— A{Saa' + l + Saa'-l)- 



(3.3) 



Here = Px/^n = Px/'m + eBvl/m. The Hamiltonian 
TYq contains three terms; a quadratic term, a linear term 
and a hopping term. If the quadratic term can be ne- 
glected, we can use as a quantum number, the eigen- 
values of Tig are zero and the eigenfunctions are expressed 
in terms of Besscl functions 



_ g-l/2^ip^z/h+i(p^-eBvl)x/n 
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2A 



zFl — eBvrl 



(3.4) 



The superscript S here means that this is a Stark-hke 
solution, and the in-plane momentum p = (pz,Px). The 

size of this state is |j/ — a| ^ A/\eBvxl — eFl\. 

We can neglect the quadratic term in Eq. ( |3.3| ) if the 
obtained state is not too extended. Comparison of the 
quadratic term with the linear term leads to the condition 
\i^-a\ < \eFl-eBvJ\/[ieBl f/m]. Therefore, the Stark- 
like wave function, Eq. (3^) is justified when 



leBvxl - eFll > hflB 



(3.5) 



where the cyclotron frequency in our anisotropic case 
is given by Eq. (1.14). In the Stark- like representation, 
the Hamiltonian Hq becomes diagonal and looks like a 
Stark ladder 



(Ep - eFlu)S,,, 



(3.6) 



where Ep — {p^ +p1)/(2m). This new interpretation of 
the stat es far above the "periodic potential" is shown in 
Fig. ^(a). 

The wave functions of the syste m a re of Landau lev- 
els type if the linear term in Eq. (3^) is much smaller 
than the quadratic term. This happens in the na rrow 
strip of the phase space Px^Pz, where condition Eq. (3.5) 
is violated. Inside this strip px is not a good quantum 
number, and it is convenient to define the continuous 
quantum number Y = mF / {eB"^) — Px/ {eB), which is a 
position of electron orbit (guiding center). Therefore we 
don't need the index v, and it will be used for labeling 
Landau type states in the rest of this section. Diagonal- 
ization of Tia leads to 



(p.,rnF/B) " sFY + HVLBiv + \) + K{Y)^6uu' , 

(3.7) 



where A^{Y) is a periodic function of Y, Aj^(F -I- /) = 
A^{Y); this function descphes the spreading of Landau 
levels into Harper bands .l3cj The cyclotron frequency 
in Eq. 1^) is defined as in Eq. ( |1.14| ). The widths of 
Harper bands can be estimated from the quasiclassical 
expression for the tunneling amplitude between neigh- 
boring Brillouin zones: 



Var A^ ^ HUb exp-^ 



4A 



i + cos(t) hnB{v + i/2) 



4A 



dt 



(3.8) 



where integration has to be performed over positive val- 
ues of the expression inside the square root. The number 
of such states is limited by the condition 



hnB{v + l/2) < 4A 



(3.9) 



V + 1 




FIG. 3.2. Approximation by Stark type states is not valid 
in hatched area only. See explanation in the text. 



and for the strong enough magnetic field there are no 
Landau type states at all. 

The width of Landau type states can be estimated from 
the usual Gaussian shape of the linear oscillator wave 
function as \a — Y/l\ ^ Av /{UVIb)- These states are 
more extended in the y direction than the Stark-type 
states and can give important contribution to the vertical 
current under very special conditions. 



3.1.2. Orthogonality conditions for Stark type wave 
functions. 

Stark like states constructed above are very important 
for further calculations and we explain now more details 
about the approximation which has been made. For fixed 
V, the Stark like state is a mixture of Wannier states from 
different wells a. The energies of these Wannier states 
are in the diagonal part of the Hamiltonian Tio minus 
pl/{2m) 



m{F/Bf \ l ( V, 



F/B 



2 \F/B 



a — V 1 f a — V 



C 



C 



(3.10) 



where C — mF/{eB'^l). We wrote this exp ression in 
dimensionless form, and plotted it in Fig. |3.2| for C = 2 
and a = i^, i/±l. So, in this figure one can see the electron 
energy in three subsequent wells. The vertical shift of the 
parabolas is due to the electric field and the horizontal 
shift is due to the magnetic field. In dimensionless units 
both shifts are 1/C. 
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Tunneling between subsequent wells mixes Wannier 
states. Two examples of such mixing are shown in 
Fig. |3.2| . Three states marked by filled circles in the ver- 
tical line v,j; = const form a new Stark-like state with 
energy close to the energy of central Wannier state. In 
order to get wave functions Eq. (3^) we hav e to n eglect 
the last term on the right hand side of Eq. (3.10). One 
may think that this approximation is valid under the con- 
dition C ^ 1, but that is not true. This appro xima tion is 
not valid for states in the hatched area in Fig. |3.2| , which 
is given by the inequality opposite to Eq. ( |3.5D . 

Let us show that Eq. (3.5) provides a condition for 
orthogonality of the wave functions Eq. (3.4). We have 



/"dr^Vfp(a,r)V'^',p'(a,r 

a 



S 



5{p'^ - Pz)S{p.^ - eBvl -p'^ + eBv'l) 

2KeBl{v'^ - v^) 
[eFl ~ eBvJ][eFl - eBv'^i 



S 



(3.11) 



We would like also to check the orthogonality condition 
by integration over momenta and summation over index 
ly. That is impossible to do, however, because of the pole 
near — F/B. 



3.1.3. Matrix elements in the Stark representation. 

The operator of the electron velocity in the y direction 
in the Wannier representation has the form 

Vy = {iAl/h){Sa^a' + l - Sa^a'-l)S{x - x')S{z - z') . 

It becomes in our Stark-like representation 
. (27r?i)2 



S 
iAd 



TP. = T 
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eFl — eBvrl 



(3.13) 



The total Hamiltonian contains Tio, and a scattering 
potential, U. We consider that the main scattering mech- 
anism is impurity scattering or scattering on the other 
type of static potential. The overlap of electron wave 
functions in different wells is small and only diagonal el- 
ements with respect to the wells, U^^,, will be taken into 

account. Matrix elements of the operators U in the Stark 
representation have the form 

TTUv' ST^ tPx tPx Tj-aa ( ^ A ] 



and Jfilc, are defined above. The product of two matrix 
elements averaged over the impurity configuration is 



X ^{Px ~ eBvl — qx + eBpl + p'^ — eBv'l — q'^ + eB/i'l) 



Y \ ^ TPx 7<?x tPx T^x 

^ 2_^'^v-a'^ p-a'^v'-a'^ p.'-a ' 



(3.15) 



where iV/ is the sheet density of the impurities and U ie 
the matrix element of the single impurity potential. 



3.1.4. Equation for density matrix 

We will assume that the Fermi energy and tempera- 
ture are much larger than the level spacing of TYq. This 
assumption allows us to derive a kinetic equation for the 
electron density matrix. We will make use of the Keldysh 
technique and derive this equation in theisame way as the 
Boltzmann equation is usually derived.L3 In the Keldysh 
technique the kinetic equation results from the Dyson 
equation, which can be written in two equivalent forms 



Gn2*Gl2 = (7z + 2130320-2 



(3.16a) 
(3.16b) 



The matrix Keldysh Green function, G12, depends on two 
sets of variables, {i/, p,i} and a sum and integral with 
respe ct to t he variables with the subscript 3 is implied in 
Eqs. ( 3.16a ) and ( 3.16b| ). The operators on the left hand 



sides, Gqj^^ and Gno^ are obtained from 



^02 



(3.16c) 



by substitution sets {vi,pi,ti} or {v2,P2,t2} cor- 
respondingly. Definitions of the Green functions 
Eqs. (3.16) contain also the Pauli matrix cr^. 

Green functions depend on two times, ti and ^2- In- 
stead of them the difference, ti — t2, and the sum, 
t = (ti + ^2)72, times can be introduced. The charac- 
teristic values of the difference time is of the order of h 
divided by the characteristic energy, i.e. the maximum of 
the Fermi energy and temperature. The time t character- 
izes much more slow variation of occu pation numb ers. So 
that in the right-hand sides of Eqs. ( 3.16a ) and ( 3.16b| ) 
all functions can be considered dependent on the same 
time t and the integration with respect to is reduced 
to the integration with respect to difference times. Then 
the Fourier transform with respect to the difference time 
leads to the Green functions depending on t and a fre- 
quency uj, which is a spectral variable. 

The explicit form of the matrix Green function as it is 
defined in Ref . ^ is 
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^ = ( G+- 



+ C G-+ 



(3.17) 



The same pattern holds for self-energies. The differ- 
ence and h alf of the su m of (— +) matrix elements of 
Eqs. (3.16a) and (3.161) has the form 



[n'^,G- 



^13 "-^32 
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"-^13^32 "-^13 ^32 ■> 



(3.18) 
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■-^13^32 



^13^^32) 



(3.19) 



For the calculation of self-energies we need also advanced 
and retarded Green functions. Equations for them are 
obtained from Eqs. (3.16a) and ( ^.16b ) 



hu;G\2--{n'^,G^}^2 



l + -{E^G'-}l2 



Ga y^r* 
12 — ^21 • 



(3.20) 

(3.21) 



Here [a, h]= ab — ha and {a, h\ = ab + ba. 

The retarded Green function can be found from 



Eqs. ( |3.20D and ( |3.24| ) (see below). Equations ( |3.19| ) and 
( p6| ) (see below) are linear and uniform with respect to 
G^jJ(p) and t his fu nction can not be found from them. 
Actually, Eqs. ( 3.19 ) and ( 3.26 ) can be used to find only 
the dependence of G ^ on uj. The integral of G ^ with 
respect to to, the density matrix 



duj 



2m 



(3.22) 



has to be determined with the help of Eq. (3.18). The 



equation for the density matrix can be obtained by the 
integration of Eq. (3.18D with respect to lo, 



duj / 
2^ V 



S'G 



-+ 



-G" + G''S-+-G"+I]" 



(3.23) 



The super lattice is uniform in the x, z plane and one 
may expect Green functions as well as self-energies to 
be diagonal with respect to p. However, this is not true 
because we are looking for a solution, which describes 
electric current, and therefore the density matrix has to 
contain terms proportional to the matrix elemen t of the 
velocity y-component. This matrix element Eq. ( ^.12| ) is 
not diagonal with respect to p, and therefore all elements 
of the Green function are non-diagonal too. 

We consider the case of weak scattering when the en- 
ergy uncertainty due to the scattering is much smaller 
then the width of the in-plane electron energy distribu- 
tion. Then self-energies can be calculated in the first 
approximation. For simplicity we assume that electrons 



in different wells are scattered by different impurities. 
This implies that the screening radius is smaller than the 
period of the superlattice. Then 



(3.24) 

E^p^,,p,(c.)=I]rp;,vH (3-25) 

(3.26) 

These self-energi es ha ve to be substituted into the right 
hand side of Eq. (|3.2S ) which can be solved together with 
Eqs. (|3l8|) and ( jslgf ). 

The next important approximation is the weak scat- 
tering limit, which allows one t o negl ect te rms contain- 
ing self-energies in Eqs. ( 3.19 ) and ( 3.20 ). Therefore 
we are going to take into account scattering in the so- 
called Born approximation. Computation of the right 
hand side of Eq. ( |3.23| ) is cumbersome and we put it into 
Appendix pT^. 



3.2. CALCULATION OF THE CURRENT 

In order to compute the current we have to find a sta- 
tionary, translational , and gauge invariant solution to the 
kinetic equation Eq. (3.23) with the right hand side given 
by Eqs. (3.47). Equation (3.23) is not gauge invariant, 
and the density matrix Ppp, also is not gauge invariant, 
because wave functions were chosen in the specific gauge. 
However, in Wannier representatiep we can introduce 
translational and gauge invariantO counterpart of the 
density matrix pa-a' (q) 



piq(r-r') 



'(q)- 

(3.27) 



Here p"^ is not an invariant density matrix defined in 
the Wannier repre senta tion. The first task is to convert 
both sides of Eq. ( 3.27 ) from Wannier to our Stark- like 
representation. We have to make another approximation 
in order to perform t his t ransformation. 

We will solve Eq. (3.23) only in the case when 



A, leBv^ -eFll-^ T^ 



(3.28) 



where Te is the effective electron temperature, which 
characterizes the width of the electron energy distribu- 
tion. (We do not assume that a real distribution is the 
equilibrium one with the temperature Tg. This quan- 
tity is used only for estimates.) In the case of an effec- 
tive energy relaxation Tg T. In the case of low tem- 
perature and an appreciable heating of the electron gas 
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Te > T. The estimate of T^ under different conditions 
can be found from the balance of the heating ~ jF and 



coohng — see Eqs. (3.39) and (3.4C) below. In the 
case of zero magnetic field the answer can be found in 
Ref. 1^. 

Under the condition Eq. ( 3.28| ) we can convert p"^ 
and pQ_Q'(q) to p^ ^, an d pi,^^'{q). Equation ( |3.27 ) is 
converted into Eq. (|3.48), see A ppendix 3.4 for more de- 



tails. Substit ution of Eq. (3.48) into the quantum kinetic 
equation Eq. (3.22) and averaging of the scattering prob- 
ability over energy surface leads to the following resulting 
equations 

hlEPoiE) + k^V{E)[R{E) + 2A^] = , (3.29) 



±i{eBvxl — eFl)p±i{p) ± iAeBvxl 



dE„ 



-r(i?p)p±i(p)-Ar(i?p) 



dpa 
dEr. 



(3.30) 



R{E)= j ^<5(£;-i?p)[pi(p) + p-i(p)] (3.31) 



where Ie is the operator of energy relaxation. Calcula- 
tions are actually performed in Appendix 3.4. Condition 
Eq. ( ^.28 ) allows us to keep the Stark- like levels differ- 
ence — z/' = 0, ±1. Quantum uncertainty of the energy 
levels due to sc atterin g T{E ^) is defined in Eq. (3.53). 
Equations ( 3. 29] ) and (3.30) are very remeniscent of the 
Boltzmann equations for the parts of the electron distri- 
bution function even and odd with respect to the electron 
momentum in the t heory of hot electrons in wide band 
semiconductorsOCS. 

The main energy relaxation mechanism in ths-^iAtrow 
band superlattices is acoustic phonon scattering.EaH We 
have 



niEP^{E) = ^Q{E) 



where 



'^Podl 



1 + 



SEmdt 



, (3.32) 



(3.33) 



Po is the crystal density, m is the in-plane effective mass, 
dw is the width of a well, and S is the deformation po- 
tential. 

Equation ( |3.30D can be written as 



P±i(p) = -A 



TjEp) ± ieBvJ dpp 
T{Ep) ± i{eBvJ - eFl) dE^ 



(3.34) 



This solution has a pole at = F/ B, which is smeared 
by scattering. The width of this smearing is F/ (eBl) and 
it has to cover the hatched area in the diagram Fig. 
The condition justifying our whole theory is, therefore. 



ngTp < 1 



(3.35) 



wher e ^ h/T{Ep). Equation ( [3.29 ) with the help of 
Eq. ( 3.32 ) is reduced to 



Q{E) 



Po{l~Po)+T 



dpQ 
dE 



where 
R{Ep) 



d<P 



2A'T{E)RiE)^ = 0, 
(3.36) 



{eFiy^ 



27r T^iEp) + {eBlpcos{(l))/m 



eFiy ■ 
(3.37) 



Integration can be performed, the same integral was cal- 
culated in Chap. ^, see Eq. (2.23). 
The solution to Eq. ( |3.36| ) is 







|exp 


\f 




Jo 



dE 



lo T + 2K^T{E)R{E)/Q(E) 



-c 




where C is a normalization constant. For zero mag- 
netic fiel|d this solution was obtained by Laikht man and 
Miller. E£l Even without an exact calculation Eq. ( |3.38 ) al- 
lows us to estimate the importance of the heating of the 
electron gas near the current maximum. We can estimate 
the momentum and energy relaxation times as Tp ^ fi/T 
and Tg ~ hTe/Q respectively. Near the current maxi- 
mum, when eFl ~ T,eBlvF,eBlvT, we have i? ~ 1, and 
the width of the electron energy distribution Te in the 
case of weak heating and the condition for weak heating 
are 



(3.39) 



Because ^ Tp, the condition Eq. ( 3.39 ) is satisfied only 
for A < T. For the band width 4A ~ 0.5 meV the last 
condition is satisfied for temperatures about 30 K and 
higher. This estimate shows that for lower temperature 
or for a wider band one can expect an appreciable heating 
of the electron gas. In such a case 



AVe 

Ttj, 



A2 



> 1 



(3.40) 



Now it is possible to justify Eq. ( 3.50| ). In the estimate 
we assume that eFl, \eFl — eBvp] ~ F because this re- 
gion of the electric field is close to the current maximum, 
and F^/A^ ^ '^e/'Tp because in the opposite case the reso- 
nance tunneling is smeared so much that the width of the 
resonance can be of the order of the separation be twee n 
the levels in a well. Then for weak heating, Eq. ( 3.39 ), 
Te ~ T > A(t^p)1/2 > A,F. In the case of strong 
heating, Eq. f^), F/T^ < (T/A)(Tp/Te)i/2 < i and 
A/Te - (T/A)(7v/7v) « (r/A)(rp/re)i/2 « l; that is, in 
both cases Eq. ( |3.2^ ) is satisfied. 
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For the uniform electron distribution, the e xpression for current density is proportional to the trace of the density 
matrix times the velocity operator, Eq. (3.12) 



2ieA f dp , , 4eA^ f dp 



(27ra)2 



T{E^)eFl 



h J {2TTh)^T^{Ep) + {eBvJ-eFl)^ dEp 



(3.41) 



This expression shows that Ohm's law is satisfied for eFl, eBvpl ^ T. The collisionless current is recovered in the 
opposite case, and 



J 



(3.42) 



This resul t is s imilar to Eq. (2.26). The important difference is the form of po; in Eq. (3.42) it has to be substituted 
from Eq. (3.3^), whereas in Eq. (2.26) we use the equilibrium distribution function /q. 
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For weak scattering, terms containing self-energies in Eqs. ( 3. IE ) and ( 3.2C ) can be neglected compared to hu and 
solutions to them are 



2 (5^^'(5(p - p') 
yp;iy'p'l^7 — T — ; — -ET, — --7: 



Glp.yp,{uo)^{2^nY 



(27r 



, V T'^ ff 



f)io- E^ + eFlfj. + iO 



2 2 

dqdq' jjv^jjf_,'u' ^ ( ^ Eg + Eg 



M + M 



(3.43) 

(3.44) 

(3.45) 
(3.46) 



All four terms on the right hand side of Eq. (3.23) can be expressed explicitly in terms of the density matrix: 



dhJ 

2^ 



(27r?i) 



\JJ [2^1^ 
Px ygx jqx jPx-eB{u-p,')l 



{Pz,Px-eB{iy-p,')l);p' 



EjPx T?x T^x TPx-e 



(3.47a) 



Z_-/Q a' —Q. 



-p'^-eB{u'-^')l jq^ jq^ jp'^ 



^ fi—a" fi—a'^ ij' —a 



Ep + S(pi,pi-ei3(.'-,,')0 ) /2 - eFZi^ -E^ + eFlp - lO 



(3.47b) 



/ (2nh)'^ f'(q,+p,,q,+p^~eB(u-t,)l)-(q,+p'^,q^+p'^-eB(y'-p.')l) 

mm' 



^ -^(g.+P.,gx+Px-e-B(»^-M)0 + E{q.+Pi,q^+p'^-eB(i^'-t,')l) _ ^ ^ _ M + ^ ^^^^z _ ^gl 
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(3.47c) 



+ E, 



{q.+p.,q^+p^-eB(l^-^i)l) ^(q^+ pj^q^+p'^-eBjiy' -f^'W _ ^ _ cFl ^ ^ + cFlv + iQ 



EtPx T«x+Px-e-B(i^-At)' jq^+p'-,~eB(v' -p.')l jp'^ 



(3.47d) 



and the expression for the commutator on the left hand side of Eq. ( |3.23 ) does not change 
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The derivative of the density matrix with respect to is proportional to l/^/mTl in the high temperature hmit. 
The expansion with respect to this parameter allow us to obtain the expression for the density matrix in the Stark 
representation. We have from Eqs. ( p. 28 ) and ( p.4[ ) 



X 6[p^-eB 



a + a 



and after shift of p^ and p'^, the above expression becomes 

^(pl,Px+e-B(i'-i^')V2);(pL,Px+e-B(i''-i')V2) ^ (27r?i)^(5(p — p') 

rP^+eB{i.^,.')l/2 jp^+eB{i.' -1^)1/2 ( ^Fn"^"' 
>^ X^Jv-a Jv'-a' Pa-a' \Pz,Px + eBl z CBl 



aa' 



= {2T:h)^S{p - p') J^-^'"'""" "'^J'l^r^r'' '"'"'U 1 - eBl 



i2nhfSip - p') 



Px+eB(iy-i/')i/2 jp,+eB{u'-u)l/2 

d 



2 2 

v — a + ly' ^ a 



AeBl 



eBpxl/m — eFl dpx 



2 dpx 



Pa-Q'(p) 



(3.48) 



where high order derivatives with respect to Px are neglected. The commutator on the left hand side of Eq. ( 3.25 ) 
becomes 



' P^' {P^.P^+'^B{u-v')l/2):(p'^,p'^+eB(u' -u)l/2) 



{eBvJ - eFl){v - iy')p„-v'{p) + AeBl{v - v')- — <^ p^_i,/+i(p) + p,.-,.'-i(p) 

OPx 



(3.49) 



The energy differences Ep — Eq in the collision operator Eq. (3.47) are of the order of T^. Therefore, we can neglect 
the terms containing in Eq. ( 3.48| ) for the density matrix when we substitute it in Eq. (3.47). We have to expand 
first collision integrals in terms of {eBvxl — eFl)/Te too. After the expansion, sums with respect to numb ers of the 
Stark levels can be calculated explicitly. Keeping only terms of the first and the second order in Eq. ( 3.23 ) we get 



i{eBvJ-eFl)l'p^{p)+iAeBlly-^^^p^+l{p) + p^^l{p)^ ^2ttNi J j^^\Vpq\^ 



S^,oPo{^) ~ P^ip) 5{Ep - Eq) 



A 



(pi(q) + P-i(q)) - Po{ol){5u.i + + p,.+i(p) + p^-i(p) 5'{Ep - Eq) 



+ 2A25,,o Po(q)-Po(p) 5"{Ep-Eq) 



(3.50) 
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We neglected the term on the cohision integral, which is different from the first one on the left hand side of Eq. ( |3.50| ) 
only by a factor. It can be considered as a renormalization of electron charge and can be neglected. Then one can 
see that if the terms of the order of A/Te are neglected at all, the equation for po is separated from equations for 
with v ^ and docs not contain fields. This is natural because po is the distribution function in a layer and without 
tunneling it does not 'know' about the electric field. Along with the tunneling in this equation it is necessary to take 
into account inelas tic sc attering, which was not considered so far. The operator of the inelastic collisions, /inpo , can 
be added into Eq.( 3.50) without making use of the Keldysh technique. 

Equations ( 3.50 ) with 7^ show that ~ (A/Tg)!"! po, i.e., p|j,| with Iv] > 2 can be neglected and we come up 
with the following equations: 



+ 



2'kNi 

n 

2-kNiK 

n 



(27r?i)2 



{2t:K) 



VpqP[po(q) - PoimSiEp - E,) + 2A^S"iEp - E,)] 

2-|VpqP[pi(q) + p-i(q) + Pi(p) + P-iip)]S'{Ep - Eg) + Lpo{p) = 



gBv I — eFl d 1 

± i P±i(p) ± AeBl — po{p) = --T{Ep)p±i{p) 



2ttNiA 



dq 



{2nhy 



\Vp^\'{po{p)-po{cL))6'{Ep-Eg) 



where for an isotropic energy spectrum and scattering (i.e. for |V^qP depending on |p — q| ) 

T{Ep) = 2nNj J -^^\Vp^\^S{Ep - E^) 



(3.51) 



(3.52) 



(3.53) 



depends only on the ene rgy. For zero electric field Eqs. ( 3.51 ) and ( 3.52| ) have isotropic solution. 

The first term in Eq. (3.51) describes an elastic relaxation in separate wells and is the largest one. If all other terms 
are neglected it leads to a distribution function depen ding o nly on the energy, i.e. po(ii)-|j^po(£'p). An equation for 
this function can be obtained by the averaging of Eq. (3.51) with respect to the energy.ESO For example 

dp 



2-Km 



5{e ~ Ep) 



2'kNi 



(^^l^pql'[pi(q) + P-i(q) + P1(P) + P-i[p)\5\Ep - Eg) 



S{e- Ep)[pi{p) + p-i(p)\ 



2'Kra 



de h 



dq 



{2'Khy 



\Vp^\^5{e~Eg 



/^^(.-i..)[pr(p) + P-r(p)]^^^ 



d 
de 
d_ 

de 



T{e)R{e) 



where R{E) was defined in Eq. (3.31). The elastic relaxation is averaged out. On the right hand side of Eq. ( 3.52| ) 
we have 



27rA^/ 

h 

2nNi 

h 

2nNi 



dq 



{2TTh)' 



\Vp^?{p,{p)~p^{q))5'{Ep~Eg) 



de 



dq 



(271^)2 



|V^pql'(po(i?p) - PQ{e))5'{Ep - e)S{Eg 



I deS{Ep - £)^ I ^^\Vp^\^{p,{Ep) - po{e))S{Eg - e) 



rjEp) dpo 

h dE„ 



Chapter 4 

Longitudinal magnetoresistance 



Classical longitudinal magnetoresistance of superlat- 
tices is calculated in the framework of a model which 
includes fluctuations of barrier conductivity. We found 
that the result depends very significantly on the fluctua- 
tions correlation length. We also found that fluctuations 
of the electron potential are not uniform along the super- 
lattice, and depend on the superlattice length. A good 
agreement between theory and experiment is obtained. 



4.1. INTRODUCTION 

In this work we consider the vertical longitudinal mag- 
netoresistance (LMR) of a superlattice; that is the mag- 
netoresistance in the geometry when both electric and 
magnetic fields are along the growth direction. Purely 
classical (i.e., without any quantum effects) LMR was 
observed many times in experiments, but only recently 
has a qualitative explanation been suggested. It is ob- 
vious that in an ideal superlattice, classical LMR has to 
be zero, because the magnetic field does not affect elec- 
tron motion parallel to it. For this reason experimen- 
tally observed LMR (Refs. |l|, H, ||, ||, ||, g ^ has 
not been explained for a rather long time;. The qualita- 
tive explanation suggested by Lee et aZ..Ej attributes this 
to nonuniform fluctuations of the superlattice barriers 
width. 

We present resistance calculations for a superlattice 
with nonuniform barriers. We consider the case of a 
narrow-miniband superlattice when the vertical trans- 
port can be considered as sequential tunneling. Each 
barrier in this case can be characterized by a conductiv- 
ity fluctuating around some average value. The opposite 
case of wide-band superlattices where an electron tun- 
nels across a few barriers between two successive scatter- 
ing events seems to be less interesting. The effect of the 
barrier width fluctuations is averaged out as a result of 
tunneling across a few barriers. 

The qualitative picture of the longitudinal magnetore- 
sistance of superlattices with nonuniform barriers sug- 
gested by Lee et alS3 is as follows. A current across each 
barrier is larger in places where the conductivity is larger. 
If high-conductivity regions of adjacent barriers are not 
positioned against each other, then nonuniform currents 



across barriers induce in-plane currents between barri- 
ers. The magnetic field perpendicular to the layers brings 
about a transverse magnetoresistance, reducing these in- 
plane currents. As a result the current across barriers 
cannot pass through places with maximal conductivity. 
In this way the magnetic field in the growth direction 
increases the superlattice resistance in this direction. 

The effective conductivity of a spatially inhomoge- 
neous medium has been considered many times iit-the lit- 
erature; see, e.g., the review paper by Landauer.llj A su- 
perlattice is just another example of such a medium with 
a specific geometrical structure of the inhomogeneities. 
We consider this problem for weak fluctuations of the 
barrier conductivity (the exact parameter will be shown 
below). We also assume that the conductivity fluctua- 
tions of different barriers are not correlated. These as- 
sumptions allow us to obtain an analytic expression for 
the superlattice resistance. 



4.2. PERTURBATION THEORY FOR 
POTENTIAL FLUCTUATIONS 

The superlattice consists of -I- 1 wells separated by 
N barriers, and the electric potential in the i^th well is 
(^i/(r ), w here r = {x,z) is the in-plane coordinate, see 



Fig. LI. The electric current ji,^^+i(r) from well v to the 
well J' -|- 1 is given by Ohm's law. 



31^,^+1 



01^+1 ) , 



(4.1) 



where cr^^+i(r) is the conductance per unit area of the 
barrier following the vth well. The formulation of the 
problem will be completed with the charge conservation 
law 



(4.2) 



where (tV(/)^ is an in-plane electric current in the vih 
well, and (j is a two-dimensional conductivity tensor of 
the well. We will assume that this tensor depends on 
the magnetic field but does not depend on coordinates. 
We will also assume that the conductivity in the wells is 
isotropic, so that Gxx = Ozz = c", and Uxz = —(^zx- This 
assumption gives 
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FIG. 4.1. Measurements of LMR in superlattice. 



(4.3) 



that is, the Hall con duct i vity does not enter into the 
problem. Equations (4.1)-(4.3) have to be solved with 
some boundary conditions. Wc will assume that poten- 
tials in the first and last wells are independent of r due to 
the presence of the highly doped uniform plane contacts, 
(j)o = NU =const and = 0. 



Equations. (O) and (42) can be solved by means of 
perturbation theory with respect to the fluctuations of 
^vv+i- the Fourier representation, 



5ay^^e 



(4.4) 



where Sa^.^i is considered a small quantity and (5cr,y,o = 0. 
The conductivity fluctuations are assumed to be uniform, 
with a correlation length much shorter than the super- 
lattice plane size, so that 



(4.5) 



where () means an ensemble average over all possible fluc- 
tuation configurations. The function is inversely pro- 
portional to the barrier area. These conductivity fluctu- 
ations induce fluctuations of the potential 



(4.6) 



The Fourier transform of Eqs. (4.1)-([0) can be lin 



earized with respect to fluctuations if q 7^ 0, 



V+l,q 



(4.7a) 



The second-order terms have to be kept in the same equa- 
tions for q = 



q 



(4.7b) 



The solution to Eq. ( [1.7a ) with the boundary condi- 
tions S(j)o^q — (507v,q = can be expressed in terms of the 
Green function, 



J_ sin(7rjV/jV) smjirju'/N) 
^yyi! ]M 2-.^ cosh(a,) - cos(7rj77V) 

1 



(4.8) 



sinh(a5) sinh(agA^) 

^ sinh(agZ/) sinh[ag(iV - ly')] , v < v' , ,^ 
sinh(agZ^') sinh[ag(A^ ~ v)] ^ v > u' , ^ ' ' 



where 



cosh(ag) = I 1 



(4.10) 



One has 

H.,^ - uY^G^.A^i)^-^^^^^^^^^^ , (4.11) 



Af-1 



with q 7^ 0. The derivation of the abov e ex pression for 



the Green function is given in Appendix 4.5. 

Physical properties of the result Eq. ( |4.li ) can be 



seen from the average value of the potential fluctuations 
squared. 



N-l 



(|^0.,qP) = f/'e,E^-^-'(9) 



2G^y{q) - Gyy-i{q) - Gyy + i{q) 



q d 



(4.12) 



where G^^^jq ) is g iven by Eq. ( [4.9[ ). 

Equation (4.12) describes the increase of the potential 



fluctuations from the contacts toward the middle of the 
superlattice. For a limited region of (?, three situations 
are conceivable. The first is the case of strong in-plane 
conductivity, when the potential fluctuations are limited 
by in-plane currents, and (|50i/,qP) ,1^ j^O,N nearly does 
not depend on v. The second is the opposite case, when 
in-plane currents are not important, and the fluctuations 
of the potential are similar to the fluctuations in a se- 
ries of random resistors. In the third intermediate case 
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the fluctuations increase with the distance from the con- 
tacts, but in the internal part of the superlattice they are 
limited by in-plane currents. 

The first case is r eali zed under the condition 
sinh(ag) > 1; then Eqs. ( |49| ) and ( |4.12| ) give 



(4.13a) 



This result is independent of v and inverse proportional 
to the in-plane conductivity. In the second case, Uq <C 
and we have 



/IXJ. |2\ TT^C ~ 
(|'50z.,q| ) = U 



(4.13b) 



In the third intermediate case, 1/A^ ^ <C 1, and we 
have to consider separately the internal region of super- 
lattice and regions near the contacts, 



aqiy <C 1 



{\5(b,,^\') = U'^q X { N~u , aq{N - i/) « 1 



otherwise 



(4.13c) 



In order to calculate the correcti on to the avera ge po - 
tential, we have to substitute Eq. (4.11) into Eq. ( 4.7b| ). 
This leads to an equation for 5 6^.0, which should be av- 
eraged with the help of Eq. (4.£). The solution to the 
obtained equation is 



K.o) = -uY, 



sinh[a,(7V - 2iy)] 2v - N 



sinh(agiV) 



N 



(4.14) 



The v dependence of the averaged potential, [N — v)U + 
{54>u,o), becomes more smooth near the contacts; that is 
near the contacts the electric field is weaker. 

One can prove that the perturbation theory developed 
here is justified if the fluctuations of the potential are 
small. 



?) « 



(4.15) 



The substitution of the results obtained above. 



Eqs. (4.13) and ( 4.14 ), to this condition leads to 

E.i:^,«„„{i,, ,4.16) 



Here go is the characteristic wave number of the function 
^g, and the quantity 1/go can be considered as a conduc- 
tivity fluctuations characteristic correlation length. 
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The total current across a barrier, 

3^a^U-5j , 



(4.17) 



is the same for all barriers, and can b e ca lculated for 
t he fi rst barrier. Substitution of Eqs. (|4.4D, (4.7), and 
( 4.11 ) in Eq. (4J.) with v — Q, averaging over all possible 
fluctuations of the barrier conductivity and summation 
over all g, give 



'Jj =a^|t/5]e,Gi4(g) + (<50i.o)| 



(4.18) 



If 

(jWq^/a^ + 4 



N-l sinh[a,(iV - 1)] 



N 



sinh(aqiV) 



(4.19) 



This equation is the main result of our paper. It describes 
the change of the current due to barrier resistance fluctu- 
ations. The sign of 5] is equal to the sign of the current 
without fluctuations; that is, fluctuations lead to an in- 
crease of the superlattice resistance. 

We are particularly interested in the application of this 
result to a calculation of the magnetoresistance of the 
superlattices. Here we consider only a weak magnetic 
field when Q.cT ^ 1, where 



eB_ 

m 



is the cyclotron frequency and r is the relaxation time. 
In this case the magnetic-field-induced change of the in- 
plane conductivity is (tII(O) — (j^\B) « fi^r^crU (0), and 
the superlattice magnetoresistance becomes 



R{B) - R{0) 

W) 



rll d 



■53 . 



(4.20) 



In transport theory, surface roughness is often approx- 
imated by a Gaussian function. Such an approximation 
immediately gives the Gaussian form for the barrier con- 
ductivity fluctuations correlation function, i.e.. 



47rS 



(4.21) 



where S is the area of the barrier, and S is the standard 
deviation of the normali zed b arrier conductivity, w hich i s 
defined generally in Eq. (4.16). Substitution of Eq. (4.21) 
into Eq. ( |4.2(]|) allows us to evaluate the superlattice mag- 
netoresistance 



R{B) ~ i?(0) 
W) 
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i3/2 



77V, jN^ < 1 , 
1117/7, 7>1, 



77V^ < 1 , 
77V2 > 1 , 



(4.22) 



(4.23) 



where 7 = crUgQ/cr^. One can see from Eq. ( |4.23D that 
the magnetoresistance disappears for both very small and 
very large go- The reasons for this, however, are different. 
The former is a case of effectively "metallic" superlattice 
wells. They are almost equipotential planes, in-plane cur- 
rents are small, and magnetoresistance is also small. In 
addition, in this case barrier conductivity fluctuations 
are averaged out and 6j itself goes to zero. The latter 
case is a case of effectively "dielectric" planes. The high- 
conductivity regions of adjacent barriers are located far 
from each other. In this case the conductances of the in- 
plane path are small and the in-plane currents are also 
small. 



4.4. DISCUSSION AND SUMMARY 

Beside the correction to the current, the surface rough- 
ness of superlattice barriers leads to a quite unexpected 
result: distribution of an electric field along the super- 
lattice appears to be nonuniform. Indeed, for the cor- 
recti on to the average potential drop across one barrier, 
Eq. ( [4.14[ ) gives 



'Ju-1,0 



2 sinh(ag) cosh[ag(2j/ — 1 



TV)] 



sinh(agA^) 



(4.24) 



In the middle of the superlattice the sign of this quantity 
is the same as that of the potential drop without surface 
rou ghne ss, C/, and near the contacts it is the opposite, see 
Fig 12. Because of surface roughness the field becomes 
stronger at the middle and weaker near the contacts. The 
size of the contact regions is about l/a<j periods. The re- 
distribution of the field along the superlattice is not a 
large effect, but it can be stronger for more pronounced 
surface roughness. The physical reason for the field re- 
distribution is that the current is trying to go across the 
least resistive regions of the barriers. Because of the lack 
of correlation of the surface roughness in different barri- 
ers, this produces an in-plane current which makes the 
overall resistance larger. Near the contacts where the in- 
plane potential redistribution is not fully developed, this 
effect is suppressed. 

The correction to the current due to surface roughness 
strongly depends on a^q^ja^ . This parameter may sig- 
nificantly vary in experiments. Its value can be estimated 
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FIG. 4.2. Contribution of the terms with different wave 
vectors to the change of the averaged field. One can see that 
the field becomes higher at the middle of superlattice and 
lower near the contacts. The scales on all the graphs are the 
same, and the dashed lines show zeros of the field change. 
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FIG. 4.3. Comparison of the experimental data of R ef. |3^ , 
(open squares), with the theoretical prediction of Eq. (4.27), 
solid line. Combinations [R{B) - R{0)]/[ER{0)nlT^] and 
(7"go/cr^ are calculated from the experimental data summa- 
rized in Table i.4 , with two fitting parameters t = 5.3 x 10~^^ 
s and go ~ 3.3 fj,m~ . 



in terms of microscopic parameters of the superlattice. 
We can estimate a-^ ~ {me'^A'^Th~'^)[l - exp{Ep/T)] 
(Ref. |26| ), where e and m are the electron charge and 
mass, respectively, A is the transition amplitude between 
adjacent wells, r is the relaxation time, Ep is the Fermi 
energy, T is the temperature, and A, H/t are assumed to 
be much less than the maximum of Ep, T. In this case 
we have 



(TSo/a± « — --^max(£;_F,T) 



2fc2 



(Iqyh 
{ArY 



(4.25) 



where / is the in-plane mean free path. For the con- 
ductivity in this expression we used the classical expres- 
sion (jll ~ ne^r/m, where n is the two-dimensional elec- 
tron concentration. This expression, as well as the phe- 
nomenological Eq. (4.2), is correct under the condition 

<C 1. The other basic equation, Eq. ([4.l| ), is justi- 
fied only under the condition of sequential tunneli ng, i .e., 
Ti/{At) > 1. That is, the right-hand side in Eq. ( [4.25| ) is 
the product of a large factor and a small factor, so that 
all cases in Eq. ( 4.23 ) are possible. In these three cases 
the temperature dependences of the magnetoresistance 
are T, T^/^, and ln(r)/r, respectively. 

The condition <C 1 means that the characteristic 
scale of the surface roughness is much larger than the 
mean free path. The theory can be easily generalized for 
the case when this condition is not satisfied. The in-plane 
conductivity in Eq. ( |4.2| ) is a response to a uniform elec- 
tric field. If the electric field is nonuniform at the scale 
of th e m ean free path, then the current conservation law, 
Eq. (4.2), holds but the in-plane conductivity cannot be 
taken from the phenomenological theory and should be 
calculated with the help of the Boltz mann equation. The 
calculation is carried out in Appendix 4.6, and the result- 
ing conductivity depends on q. The only modification in 
the previous theory is that the expression crllg^ is deter- 



mined now by Eq. (4.42). The estimate of the magne- 
toresistance, Eq. (4.23j7which was done before for the 
phenomenological case, is replaced in the limit qol ^ 1 
by 



R{B) - R{Q) _ , 2 
i?(0) -^"-^j 



In 



9o« 



A^T- 



(4.26) 



In the calculation of the conductivity we neglect quan- 
tum corrections. That is justified when the magnetic 
quantization is smeared by scattering, flcT <C 1, or at 
high enough temperature, when hflc ^ T. 

The comparison of the our result with available experi- 
mental data is difficult because not all parameters neces- 
sary for theoretical calculations are known. Here we com- 
pare our results with the measurements of Ref. taking 
the relaxation time and the characteristic length of in- 
terface roughness as adjustable parameters. The widths 
of wells (dm) and barriers (c?_b) and measured magne- 
toresistance and barrier conductances in four measured 



samples are summarized in Table 4.4. We assume that 
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TABLE 4.1. Summary of experimental data, which was 
used in Fig. i.S . 



(A/A) 



AR/{RB') 



= N/{RS) 
S/cm^ 



(50/50) 
(20/80) 
(80/20) 
(20/40) 



0.36 
0.09 
0.25 
0.19 



0.027 
0.032 
0.042 
0.047 



1200 
5500 
4200 
14000 



fluctuations of the transition amplitude between adja- 
cent weUs, A A, result from the fluctuations of the width 
of the barrier by ±1 monolayer. The known geometry of 
the structures allowed us to calculate A, and the fluctu- 
ations of the barrier conductance were evaluated accord- 
ing to S « (2AA/A)^. The conductivities in wells were 
calculated according to (jH = ne^r/m, where the two- 



dimensional electron concentration n 



10 



17 



The dependence of the magnetoresistance on parameters 
of the samples can be written in the form 



RiB) - i?(0) 



i?(0)B2 e2 



(4.27) 



For the function in the right-hand side Eq. (4.22) gives 



/(7) 



dxe 



{jx + 4)3 



1/2 



(4.28) 



In Fig. 4.3 we show the theoretical curve for this function 
and the experimental results for four samples. The best 
fit is obtained for t — 5.3 x 10"^^ s and qo = 3.3 /xm~^. 
We have to note an obvious qualitative and good quan- 
titative agreement between the theory and experiment. 
The discrepancy (no more than 10%) can be attributed 
to slightly different relaxation times r and surface rough- 
ness characteristics in different samples. 

The results of the fitting give qqI « 0.2, which justi- 
fies the phenomenological expression for crU used in the 
calculations. The estimate of the surface roughness re- 
laxation time for a quantum well with d^j = 50 A and 
qo = 3.3 /xm~^ gives a value of 2 x 10^^° s. That is, 
the dominant scattering mechanism is probably impu- 
rity scattering, which explains approximately equal re- 
laxation times in samples with different values of d^. 
The surface roughness correlation length of 3000 A seems 
large, but even values larger by an order of magnitude 
have been reported. 

In summary, we calculated the correction to the super- 
lattice resistance due to nonuniform fluctuations of the 
conductivity of each of the barriers. Our results explain 
the classical longitudinal magnetoresistance of superlat- 
tices. We found that the magnetoresistance has a non- 
trivial dependence on the characteristic length scale of 
fluctuations; it goes to zero for both very large-scale and 



very-small scale fluctuations. The fluctuations of the bar- 
rier conductivity lead also to a nonuniform distribution 
of the electric fleld along the superlattice. The results 
of the theory give a good quantitative agreement with 
experimental data. 



4.5. DERIVATION OF THE GREEN FUNCTION 

Trigonometric sums can be calculated with the help of 
the equation 

2 cosh{X)u{v) - u{iy+l)-u{iy-l)^0 , < v < N . 

(4.29) 



The general solution to Eq. ( 4.29 ) has the form 

u{n) = Cie^" + Cae"^" , (4.30) 



and the eigenvalues and the eigenfunctions of Eq. ( 4.29| ) 
with the boundary conditions uq = = are 



cosh(Aj) = cos{tt j/N) , l<j<N~l 
"j (^) = ( 77 ) sin(7rjV/iV) . 



(4.31) 



Let us consider the function Gu^^i satisfying the equa- 
tion 

2 cosh{\)G{u, u') - G,+i,,. - G,_i.,. = (4.32) 

with the boundary conditions Go,i^' — Gisiy — 0. 
This can be expressed in terms of the eigenfunctions 
Eq. (pl|) , 



Af-1 



G 



sm{-Kju/N) sm{TTjv' /N) 
cosh(A) — cosh(Aj) 



(4.33) 



and this is the first line in Eq. (4.8), with A = a^. 

The same function can be calculated i n ano ther way 
with the help of the general solution Eq. (4.30). Appar- 
ently, Eq. (4.32) with the boundary conditions is satisfied 
by the function 



Gy,,y' = Ai sinh(Ai/) , v < v' , 
G^y = A2 sinh[A(iV - i^)] , v > v' 



(4.34) 



The only problem is with the values v = v' . The values 
of A, i?, and G 1,1, 1,1 ca n be obtained by substitution of 
the solutions Eq. ( |4.34| ) into Eq. ( |4.32| ) iox v ^ v' ,v' ±\, 
that leads to the system of equations 

Ax sinh(Ai'') - G^^y = , 
Ax sinh[A(j/' - 1)] - 2cosh(A)G^',^' 

-^^2 sinh[A(iV -v' - 1)] -1 , 
G^,y - A2 sinh[A(iV - v')] = , 

(4.35) 
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which has a solution 



A^ = 



sinh[A(iV - v')\ 
sinh(A) sinh(AA^) ' 
sinh(A!/') sinh[A(A^ - v') 

sinh(A) sinh(AiV) 

sinh(Ai^') 



sinh(A) sinh(A7V) ' 
and the final expression for the Green function is 



(4.36) 



sinh \v sinh A(A^ — v') 



V <v' 



smli A sinh AA^ ' 

(4.37) 

This result is also an explicit expression for the sum in 



sinh Ai/'2 sinh A(jV - v) 
sinh A sinhAA' 



Eq. (4.33) and this is the second expression in Eq. (4.9) 



1 - Jn I 2— sin 



(4.42) 



For the case ql <C max(ilcT, 1), this equation gives the 
usual classical result, 



1 + f72T2 



(4.43) 



4.6. WAVE- VECTOR-DEPENDENT IN-PLANE 
MAGNETOCONDUCTIVITY 



The Fourier transformation of the Boltzmann kinetic 
equation for the two-dimensional electron gas in the con- 
ducting layer is 



Of 



iq ■ v/q + e[v X B]^ - zeq • v0q^ = , (4.38) 



dp 



where f{E) is the equilibrium distribution function, /q 
is the Fourier transform of the distribution function per- 
turbation, 0q is the fluctuation of the electron potential, 
and B is the magnetic field, which is c onsidered to be 
perpendicular to the layers, see Fig. 4.1 i.e., q _L B. By 
introducing q = (g, 0,0), v = (w cos(6'), 0, u sin(6')), and 
B = (0, mfJc/e, 0) we reduce the kinetic equation to 



de 



— h iqv cos{9) 

T 



-iqve4>ci cos(e')^ 



(4.39) 



It is easy to solve this first-order differential equation, 
and the answer is a function of v and 9: 



/q 



e0q 



dE\ 



10 

. qv 



dO' 



X exp( — + [sm{d) ~ sin(6l' + d)] 



(4.40) 



For Eq. (4.2) only the current divergence is necessary, 
2dp 



«q - Jq = 



(27ra)- 



; ieq- v/q = crllq^c; 



(4.41) 



The last equality is the definition of crU. An integration 
with respect to angles in Eq. (|4.41 |) can be performed, 
but the integration from Eq. ([4.40[) still remains 



Chapter 5 

Theory of high-field-domain structures 
in superlattices. 



A number of experimental works provide evidence for 
the existence of the high-field domains in superlattices 
when the applied voltage exceeds some critical value. A 
theoretical description of the structure of such a domain 
is developed. We confine ourselves to the case of narrow- 
band superlattices, where electrons are strongly localized 
in the wells. We find that the minimum length of the 
high-field domain can be larger than one superlattice pe- 
riod. The maximum current in the oscillating part of 
the I-V characteristic can be significantly smaller than 
the value of the current at the voltage where the first 
instability comes about. The oscillation period can be 
considerably smaller than the value corresponding to the 
energy separation between the first and the second level 
in a well. For the case of the domain formation at some 
distance from the anode, we study the field distribution 
in the low-field region downstream of the domain. 



5.1. INTRODUCTION 

During the past 20 years a number of interesting ex- 
perimental works have been performed in order to inves- 
tigate traasaort nrimeijtigs.Qlj si.iuet]attices in the growth 
direction.QO'EllHSa'LllH'Q'EjS Under a weak ap- 
plied bias the superlattice looks like a homogeneous 
medium and exhibits Ohm's law. Near some critical field 
-Fth an instability appears and destroys the homogeneous 
state. As a result of the instability the superlattice breaks 
down into three regions: the low-field region with trans- 
port in the first mini-band, the high-field domain and the 
low-field region where electrons are injected into the sec- 
ond miniband from the high-field domain and then relax 

An electron 



down to the first miniband; see Fig. 5.1 



can nMyfi-LOpOA in the second mini-band before it drops 
down,E3EZlE3 because the inter-subband relaxation rate 
is relatively small. A further increase of the applied bias 
leads to an expansion of the high-field region and the 
current exhibits an oscillatory behavior. The period of 
this oscillation can be associated-Awith the intersubband 
space but generally it is smallerE3'EJ. Under higher bi- 




1^ = I'll 



FIG. 5.1. Regions of different conductivity in the superlat- 
tice. I and III are low-field domains, II is a high-field one. 
Dashed lines show the level positions. Levels are broadened 
due to scattering. The second levels in region III form a mini- 
band and the long dashed lines show its edges. Arrows show 
the hopping of electrons between the levels. In region III most 
of the electrons move in the second mini-band. 



ases upper minibands become involved in the transport 
process. 

There is no general physical law forbidding domain for- 
mation at any place in the superlattice. However, in the 
undoDfid superlattices a domain appears naturally at the 
anodcES, see Figs. 1.3 and 5.2(a). In the doped superlat- 
tices domain can be formed in the middle of the super- 
lattice. Fig. 5.1, or near the cathode. Fig. 5.2(b). In the 



last case one has to observe a significant increase of the 
current after the instability point. 

A phenomenological model that described a superlat- 
tice by an jequivalent electric circuit was suggested by 
Laikhtman.E3 In this model each barrier was replaced 
by a nonlinear resistor parallel to the capacitor. The 
model explained current oscillations and the hysteresis 
usually observed in the experiment. Prengel, Wacker and 
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also discussed. 

The physical processes characterizing transport in the 
supe rlat tice are described briefly in the next section. In 
Sec. 5.3 we derive the equation for the hopping current 
between two levels in the different wells. In Sees. 5.4, 5.E, 



V = VII 



FIG. 5.2. Particular cases of the potential distribution from 
Fig. 5.1. Domain growth from anode, (a), this is the most of- 
ten observed experimentally case. Domain growth from cath- 
ode, (b), this situation has never been reported in the litera- 
ture, however, it has very specific fingerprint on the I-V curve 
and it can be observed. 



SchollEl considered a model for a realistic superlattice 
which included electron tunneling between different levels 
in the adjacent wells and the relaxation processes inside 
one well. They obtained multistability of the current- 
voltage characteristic and various hysteretic transitions 
which arose upon sweeping the applied voltage and which 
they associated with changes in the domain size. 

The purpose of the present work is not to simulate the 
I-V characteristics in a specific superlattice but to under- 
stand the general structure of the high-field domain. A 
diffusion current induced by a charge accumulation at the 
domain boundary appears to be very important. We cal- 
culate the field and the carrier distribution in the steady 
state and get the main features of the current-voltage 
characteristic of the superlattice in some interval of the 
applied bias. The size and the position of the domain is 



and 5^ we calculate the field distribution in regions I, II, 
and III respectively. We discuss the results and make 
some comparisons with available experimental data in 
Sec. O. 



5.2. PHYSICAL PICTURE 

Let us see what happens when the bias applied to 
the superlattice increases and goes beyond the instability 
threshold. If as a result of the instability developruent a 
short high-field domain is spontaneously generatedEl then 
the field in the low-field regions is reduced by F^i/N, 
where Fa is the field in the high-field region and N is 
the number of the superlattice periods. The current in 
the low-field region is j « 2joFFth^/ {F"^ -\- i^^)' where 
i^th is the threshold field r«uid jo is the current just be- 
fore the instability pointfis. So the generation of the 
short domain causes the reduction of the current by 
6] ~ jo{Fu/FthN)^. Usually Fh is only 10 or 15 times 
larger than i^th, but the number of periods can range from 
50 to 100, and therefore the reduction of the current is 
small, Sj <C j. 

The formation of the high-field domain is accompanied 
by the accumulation of electrons in the well just upstream 
of the domain and with the depletion of electrons in the 
well just downstream of the domain. The accumulation 
of electrons in one well gives rise to a diffusion current 
upstream of this well in the direction opposite to the total 
current. Since the total current is the same across all of 
the barriers the diffusion current across one barrier has 
to be compensated for with the conduction part of the 
current. This compensation may be impossible, because 
the electric current is too close to its maximum value jo- 
That is, for a such a value of the total current a steady 
state does not exist. A steady state can come about only 
for a domain extended enough when the total current 
isn't too close to its maximum value jo. Therefore there 
exist a minimal length of the high-field domain and an 
upper limit of the total current in a steady state j*. 

Generally, after the formation of the high-field domain 
with an increase of the applied bias, the total current 
drops below j*. A further increase of the bias leads the 
growth of the current and when it reaches the j* the 
high-field domain expands by one period and the current 
drops again; see Fig. 1.2. 



The change of the potential drop across the high-field 
domain when it expands by one period is usually associ- 
ated with the energy sp|acin^|£L,|Jj|ipt|Wpen the first and the 
second levels in a wellJaOEaM'Btje It is assumed that 
in the high-field domain the first level in one well is in 
resonance with the second level in the neighboring well. 
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The number of electrons in the first n^-*^^ and second n^^^ 
levels can be found from the simple balance equation, 



.(1) 



.(2) 



,(2) 



T21 



(5.1) 



where n'-^-' + 71^2) — is the total concentration. Here 
hT /AI2 is the transition time between adjacent 



Tt 



wellsEl, r is the width of the level, A12 is the overlap be- 
tween the wave functions of the first level in one well 
and the second level in the adjacent well, and T21 is 
the relaxation time from the second level to the first 
one in the same well. If Tt <^ T21, then the current is 
en[l — exp{—Eg/T)]/2T2i. This current must be smaller 
than j* , which is not always the case. If Tt 3> T21 the 
current is ef7,tanh(i?g/r)/rt. However this quantity is 
even larger than jo = enA^-^/ max(i?F, T), where An is 
the overlap betweeji the wave functions of the first levels 
in adjacent wells .c3 

We see that under resonance conditions the current 
in the high-field domain sometimes appears to be larger 
than the maximum possible current in region I and such a 
regime cannot exist. Due to the limitation of the current 
in region I, resonance in region II is not reached. The 
current in this region is smaller than its resonance value 
for two reasons. The tunneling probability is reduced 
because of a lack of resonance and not all electrons in 
the second level in one well have enough energy to move 
to the first level in the neighboring well. 

If the resonance between the first and the second level 
in adjacent wells does not exist, the expansion of the 
high-field domain by one period requires a voltage in- 
crease smaller than that corresponding to Eg. This is 
the explanation of a small period of the jnprpntjpscilla- 
tions sometimes observed in experimentsoEiJifj'LJ. 

The well at the boundary between regions II and HI is 
depleted. The reduction of the electron concentration in 
this well corresponds to a field discontinuity between the 
high-field domain and region III. If the necessary reduc- 
tion is larger than the average electron concentration in 
a well, then the domain is located near the anode, where 
a depletion layer is formecEl. 

In doped superlattices electrons come from the high- 
field domain to the second miniband in region HI and 
relax there down to the first miniband. The relaxation 
length depends on the relation between the mobilities in 
the first and the second minibands and an intersubband 
relaxation time. The redistributing of electrons between 
two minibands can result in a field inhomogeneity in re- 
gion III. 



5.3. ELECTRIC CURRENT BETWEEN 
ADJACENT WELLS. 

In this work we consider the case of elastic scattering so 
strong that an electron is scattered in a well before tun- 
neling to the next well, at least in the first miniband. So 



for the calculation of the current we need the transition 
probabilities between adjacent wells. Since the widths of 
levels are typically much smaller than the energy sepa- 
rations between them, only the tunneling between those 
levels that are close to resonance is important. 

The general form of the transition probability between 
such levels is 



W (X 



A2 



(5.2) 



where A is the energy separation between levels. In the 
low-field region this equation describes the transition be- 
tween lowest levels and it is justified for An <C F; see 
Ref. For the transition from the first level to the 
second level in the adjacent wclL such an equation wa s 
derived by Kazarinov and Suris.H In this case Eq. (5.2) 
is justified for an arbitrary relation between A and F. 

The overlap integral A and the level width due to an 
elastic scattering F are different for different pairs of lev- 
els. The overlap integral increases with the level number 
because the penetration length of the wave function un- 
der the barrier increases with the energy. The parameter 
F in Eq. (5_^) is different for the transition from the first 
level to the first level and for the transition from the 
first level to the second level, because in the latter case 
the presence of the first level gives more possibilities for 
momentum relaxation. 

In a part of region HI electrons can travel in the second 
subband, which can be wide. Therefore in this region 



instead of Eq. (5.2) we use Ohm's law. 

The transition probability Eq. (5.2) gives the following 
expression for the electric current from the i-th level in 
the i^-th well to the i'-th level in the v + 1-th well: 
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h 

>0, 



2dp 



2FA2 



{2nhf F2 + At 



■{piEp)-p'iEp 



when the level in the i^-th well is higher and 



A 



u.ii' 
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~ h 
<0, 



2dp 2FA2, 
(27r?i)2 r2 + Al 



■{p{Ep - A^^ii>) 



-A^,u')) , 
(5.3a) 

(5.3b) 



when the level in the v-th well is lower. In these equa- 
tions Aiii is the overlap integral between electron wave 
functions of the levels i and i' in the i^-th and u + 1-th 
wells correspondingly. The energy space between these 
two levels is denoted by A^.^/. The diagonal elements 
of the electron density matrix related to these two levels 
p and p' can be considered as a function of the energy 
Ep = p'^/2m, where m is the effective mass of the elec- 
tron, since the in-plane motion of electrons is isotropic 
and this density matrix element is independent of the di- 
rection of p. Usually F is a smooth function of the energy 
and we assume it to be a constant. 

We see from Eq. (5.3) that the value of the current de- 
pends on the shape of the electron distribution function. 



5.4. LOW-FIELD REGION UPSTREAM OF THE DOMAIN. 
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Equation (5.3) is simplified in three cases. The first is 
the case of weak electron heating when the electron dis- 
tribution function is close to the equilibrium one. 

In the second case the electron gas is degenerate and 
A smaller than the Fermi energy. Then the difference of 
the distribution function in the integrands of Eq. (5.3) is 
proportional to A and the tail of the distribution function 
above the Fermi energy does not play any role. 

In the third case the electron gas is heated significantly 
so that the electron-electron scattering is very effective 
and leads to a fast relaxation of the electron distribution 
function to the Fermi function with an effective temper- 

(i) 

ature, T, and a chemical potential Q , where u is the 
index of the well and i is the index of the level. We 
should note, however, that even a strong deviation of the 
electron distribution from the Fermi function does not 
change qualitative results of the present work. 

Under these two assumptions the integration in 
Eq. (5.3) results in 



2eA? 



(5.4a) 



when the level in the v-th well is higher, and 
2eAl, r 



^ r2-fA2 



(5.4b) 



when the level in the i^-th well is lower. In these two 
equations 



(5.5) 



is the concentration of electrons in the h'-th well at the 
i-th level and go = m/Trh^. We will omit the subscripts 
of A when it does not lead to a confusion. The barrier 
for the second level is lower than for the first level, so one 
can expect that An < A12 < A22. 

The difference in the square brackets in Eq. (5.4) can 
be simplified. 



n(C<^-A,,,,) 



i^A,,,,,( — j 



(5.6) 



in any of two cases, C-A>TorC<T, A<r. 

In the case when the expansion Eq. (5.6) is used one 
can distinguish between the diffusion current and the 
conduction current. The former is proportional to the 
concentration difference and the later is proportional to 
A. Note that 



dn 



9 = 30 (1 - e" 



\/goT 



(5.7) 



is not a constant, but depends on the electron concentra- 
tion, which can be different for different levels. 



5.4. LOW-FIELD REGION UPSTREAM OF THE 
DOMAIN. 

We assume that in the low-field region upstream of the 
domain there are electrons in the first miniband only. 
Motion of electrons in the narrow miniband A <C F 
can be. described in terms of hopping between adjacent 
wells Ej. S o, a current via each barrier can be found from 
Eq. (5.4a), where i = i' = 1, and these indices will be 
omitted throughout this section. 

The electric field in this low field region is inhomo- 
geneous only near the boundary with the high-field do- 
main. The field distribution near this boundary can be 
calculated from the Poisson equation together with the 
condition that the current is the same throu gh al l bar- 
riers in this region. One can see from Eq. (5.4a) that 
the current through a barrier is a nonlinear function of 
the electron concentrations near this barrier and of the 
field in this barrier. A/ (el). For simplicity we consider 
only the degenerate electron gas. In this case dn/dC in 
Eq. (5^) is a constant; it is equal to go, see Eq. (p/^). 
This restriction is not very strong since additional elec- 
trons come to this region from region III if it exists or 
from the anode contact. Therefore the electron gas in 
region I near the boundary of the high-field domain is 
typically degenerate. 

In this case the condition that allows us to use the ex- 
pansion Eq. (5.6), is A^, < n^+i/go. Since the current 



in the superlattice with the high-field domain is smaller 
than jo, the potential drop per period far from the do- 
main boundary in region I is small, i.e. A < F. On 
the other hand, the theory is limited by the condition 
F ^ Ep, and we see that far from the domain A < Ep. 
The field in the barriers increases with the approaching 
the high-field domain boundary, however, the electron 
concentration also increases and the necessary condition 
is usually fulfilled. It makes sense to note that the neces- 
sary condition (for A^) contains the concentration down- 
stream of the v-ih barrier where it is larger than that 
which is upstream [rii, < nj/+i). 

We introduce here two quantities which can be mea- 
sured in practice. The first of them is the linear conduc- 
tivity in the low subband a. One can find from Eqs. (5.4)- 
(|j1) that a = 2e'^glkl^/nT, see also Ref. 1^. The second 
one is the critical field Fth. This field corresponds to the 
instability of the homogeneous steady state. In Ref. ^ 
it was shown that f th /el. The substitutions of these 



quantities into Eq. (5.4a) gives 



)/e9 



I l + [{^u-^u+l)lFtl,l] 



(5.8a) 



where v < ui, and vi is the number of the well between 



the regions I and II, see Fig. |5.1 . 

Here ipi, is the diagonal matrix element of the elec- 
tric potential in the z/-th well. The term proportional 
to the concentration difference on the right hand side of 
Eq. (5.8a) is a diffusion current. In region I the electron 
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concentration grows in the vicinity of the high-field do- 
main and therefore the direction of the diffusion current 
is opposite to the direction of the total current. In terms 
of these potentials Aiy^ai = e(f>i, — e(f),y^i. 

The given definition of the potentials allows us to avoid 
taking in to consideration the well polarization. This ef- 
fect is taken into account in Ref. where the integrated 
Poisson equation was derived, ft provides necessary con- 
nection between potentials and concentrations 



A4(b] + —A4n] 



{n„ — n) . 



(5.8b) 



Here C^s and e are constants, which can be calculated for 
a given superlattice, A^[/]=/(iy -I- 1) -I- f{y — 1) — 2/(i/), 
and n is the average electron concen tratio n. The second 
term on the left hand|-side of Eq. ( ^.8b| ) describes the 
capacitance of one wellc3. 

The system of Eqs. (5.8) has two boundary conditions. 
First, when v goes to — oo the difference (p^, — (j^v+i goes 
to Fool. Second, Eq. ( 5.8b ) for v 
fcrcnce 6,,, 



vi contains the dif- 
•vi-\-\ which is defined by the field in the 
high- field region, i^jj, and therefore it i s abo ut FhL Far 
from the boundary with region II Eq. (5.8a) becomes 



gF^ 



J = 



1 + {Foo/Fthf 



(5.9a) 



This equation shows that the current in region I is limited 
from above and reaches a maximum at Foo = -Fth 



i < ja 



1 



-aF, 



th 



(5.9b) 



It is convenient to introduce a dimensionless field fi, 
and displacement as following, 



HFl-Fl) 



K - K-i = 3^ _ p2 ) ("^ - ' 

where v < vj. Equations (5.8) take the form: 
- Ar'A,4/i] - (3/2)(/„)2 = , 



where 



A 



4^eV Fl-Fl 



^th ^ ^oo 



Aniens 



The boundary conditions for Eqs. (5.11) are 

lim = , 

iy—>- — oo 

2(Fh - Foo)Foo 



(5.10a) 
(5.10b) 



(5.11a) 
(5.11b) 



(5.12a) 
(5.12b) 

(5.13a) 
(5.13b) 



The coefhcient A2 depends only on the superlattice pa- 
rameters. From the definitions of the e and CeS in Ref. |26| 
one can easy get that A2 > 2. The parameter Ai depends 
on i^no and therefore o n the current through the superlat- 



tice \l cx \/l — P/ji, see Eqs. (5.9, 5.12a ). Therefore Ai 



goes to zero when the current approaches its maximum 
value. 

The analytic solution to nonlinear difference equations 
Eqs. (5.11) can be obtained in some limiting cases, but 
we consider here only one important example. Later we 
give a numeric solution in the general case. 

For Ai <^ 1,A2 the variation of fi, and hi, from well 
to well is small. Hence, the second difference ] can 
be replaced with the s econd derivative (P /dv^ and it can 
be neglected in Eq. (5.11b). The resulting differential 
equation has the solution 



1 



cosh (Aiz^/2 -f const) 



(5.14) 



The constant in Eq. 
boundary condition Eq. 




can be found from the 
The important prop- 
erty of this solution is that it is limited from above. Such 
a limitation is not connected with a small value of Ai 
but is a general property of Eq. (5.11). This limitation 
ultimately results from the limitation of the current in 
the first miniband; see e.g., Eq. (5.9a). In general, the 
upper limit of the solution to Eq. (5.11), which satisfies 
boundary condition Eq. ( ^.13a ), depends on Ai and A2. 



<T(Ai,A2) 



(5.15) 



Equations (5.11) can be reduced to the recurrent rela- 
tion 



fu-i = QUu) , 



(5.16) 



where the function Q(x) depends on the parameters Ai 
and A2 and does not depend on v. T his fu nction has 
to satisfy the boundary condition Eq. (5.13a) i.e. Q{x) 
vanishes when x goes to zero. For 2; ^ 1 function Q{x) 
can be calculated explicitly. 

Typical plots for Q{x) are shown in Fig. ^.3[ One can 
get the sequences of values of fi, by iterations of the func- 
tion G{ x). Fo r example, fiyj-2 = QiQiUi)), fui is given 
by Eq. ( 5.13f| ). These iterations are shown in Fig. 5.3 by 
dashed lines. One can see that 



T 



(5.17) 



where Q^^{x) is the function inverse to G{x), 
G{Q^^{x)) — x. The position of this maximum is also 
marked in Fig. 5.3 by a square. In Fig. 5.4 this maxi- 
mum T is plotted as a function of A^ for different values 
ofAi. 

A calculatioECj shows that \\ > 4:{dB£cs/dw^ + l)j 
where ds and dw are the widths of the barrier and the 
well respectively, e is dielectric constant in the barrier. 
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FIG. 5.3. Q{x) for diflerent values of \\ and = 20. The 
quantities ju, which are proportional to the electric field in 
the barrier i/, can be obtained by iteration of the function Q(x) 
and the dashed lines show the example of such an iteration. 
The iteration starts from the value of / that corresponds to 
the field in the high-field domain. Squares show the maximum 
value of that /. 
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FIG. 5.4. T as a function of Af is plotted for different values 



and Eeff is the effective dielectric constant in the well. 
Although Eeff is proportional to the number of electrons 
in the well, usually eoff ~ e and de dyj. Thus a typical 
value of A2 is 8-10. It becomes large only in the limit of 
extremely narrow wells or highly doped wells. The other 



parameter, can be estimated by means of Eq. ( p.l2a| ) 



It gives < e^gol — Even if the superlattice 

period is 500A A^ < 20. We plot T for A^ that ranges 
from to 60; see Fig. |5.4 The estimate for typical values 
of A^ and A2 shows th at us ually T « 1. 

The condition Eq. (5.15) can be rewritt en in terms of 
the current j by making use of Eq. (5.13b) and Eq. (5.9) 



3 < J 



where j* can be found from following equation 



r 



jo + (3T-l)Vj1-(j*)2 



-Fth 



(5.18) 



(5.19) 



The quantity j* is the upper limit of the current in 
the superlattice with the high-field domain. The c urren t 
through the superlattice is also limit ed by Eq. ( 5.9l| ), 
but j* < jo, and therefore Eq. (5.18) is a stronger re- 
striction. This restriction comes from the properties of 
the boundary between regions I and II. Indeed, the ex- 
cess of electrons in the well at this boundary caused by 
a large field gradient generates a diffusion current oppo- 
site to the current through the superlattice. Usually the 
diffusion backfiow of electrons is compensated by a local 
increase of the field. This compensation is possible only 
if the current through the superlattice is smaller than the 
maximum current in region I. 

The condition Eq. ( 5.18 ) also implies that there exists 
an upper limit of the field in the low-field region. When 
this field exceeds that limit value, the system becomes 
unstable. The development of this instability leads to 
the expansion of the high-field region, and the field in 
the low-field region decreases abruptly. 



5.5. DESCRIPTION OF THE HIGH-FIELD 
REGION. 

In the high-field region the main contribution to the 
current is the electron tunneling from the first level of 
one well to the second level of the neighboring well fol- 
lowed by the relaxation from the second level to the first 
one. The main mechanism of the energy relaxation is 
the emission of optical phonons, if the inter-subband en- 
ergy space Eg is larger than the phonon energy hflho ■ In 
this case the relaxation time T21 ranges from 0.5 x 10^^^ 
to 1 .?<^{: . . drnf ^"^'"5 on the inter-subband energy 



spacetj'0EiHl'BESEa. This is larger than the relaxation 
time in the bulk material, because the scattering prob- 
ability is inversely proportional to the transfer momen- 
tum squared and the transfer momentum in the inter- 
subband relaxation, Im^fiEg — fiQ^i^o), is larger than 
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that in the bulk material 2m[y/lE) — ^/{E — hQi^o)] ~ 
see Ref. In the case of Eg < hili^o 
the main relaxation mechanism is the electrjon-electron 
interaction and T21 is about 10^^" - 10~^ sec.EjEjH The 
transition time Tt can vary in a wide range depending 
on the superlattice parameters. This time has been mea- 
sured in 40 A/40 A and in 30A/30A GaAs-GaAlAs super- 
lattices, see Ref. and appeared about 3.6 x 10~^^ sec 
and 5.3 x 10~^^ sec respectively. In a 123A/2lA super- 
lattice Tt has been found about 6 x 10~^^ sec; see Ref. 

In the case when the first and the second levels in adja- 
cent wells are not in re sona nce one can not use the simple 
balance equation Eg. (|5.l| ). The current J12 in region II 
is described by Eq. (5.4h|) with i — 1 and i' — 2. The 



energy space between these levels A = eF^l — Eg < 0. 
The current is equal to the number of electrons that relax 
from the second level to the first level per unit time 



(2) 



.712 



= e- 



T21 



(5.20) 

For simplicity we neglect in these calcula- 
tions exp(— i?g/T) compared to unity, because usually 
Eg > T. So, the generalized balance equations become 



„(^(i)+A)_n(^) _ ^ 

n T21 



(5.21a) 
(5.21b) 



In this equation, the inverse transition time l/rt = 
2Al2T/h{A'^ + r^). 

The elimination of rt^^^ from Eqs. (5.21) leads to 



e 2Af2rn(C(i) + A) 



^ a A2 



r2 + 2A22rT2i/?i 



(5.22) 



This equation together with Eq. (5.21) describes the 
current-voltage characteristics of the high-field domain, 
i.e. the dependence of the current on the potential drop 
per period, Fl. In general, this quantity is different from 
the value corresponding to the resonance between the 
first level in one well and the second level in the neighbor- 
ing well, Eg/e. Usually it is assumed tkatitiie deviation 
from the resonance is negligibly smalE'EaES. Actually 
the levels can be considered in resonance only if |A| < T. 
However, this is not al way s the case. 

The results of Sec. show that the current in the 
superlattice with the high- field domain is smaller than j*. 
If T21 is not very large, then 



j* < en/2r2i 



(5.23) 



In the other possible case, j* « en/2T2i, the levels have to 
be in the resonance and therefore A 0. The condition. 



Eq. (5.23), means that n'^) <^ n, and the last inequality 

i(C(^) - |A|) < n or Tt > T21. 



is satisfied in two cases, 
In the first case 



\A\>T,Ef 



(5.24a) 



i. e. the deviation from the resonance is rather large. In 
the second case A^ -f E^ > 2A^2rT2i/?i. Usually F --3-5 
meV, T21 > 0.5x 10~^2 sec, A12 ^3-5 meV, and therefore 
F < Al2T2i/fi. That is, in this case 



A| » JAl2rT2i/h > F 



(5.24b) 



The se in equalities for A show that under the condition 
Eq. (5.23), the increase of the applied bias necessary for 
the extension of the high-field domain by one period can 
be considerably smaller than the resonance value Eg/e. 

The limitation of the current by the value of j* orig- 
inat cs fr om the boundary between regions I and II, see 
Sec. 5.4, When, with the increase of the bias, the high- 
field domain extends over the entire superlattice, this 
boundary disappears. Then the current jumps up sharply 
and reaches the value defined by conditions of the reso- 
nant tunneling. 



Jres 



2A2 

ZYVj2 



Tt + 2t2 



AAl2T2i/n 



(5.25) 



5.6. 



LOW-FIELD REGION DOWNSTREAM OF 
THE DOMAIN. 



In this region, if it exists, electrons are injected into 
the second miniband from the high-field domain. They 
can move-jlOOOA before they drop down to the first 
minibandEll. The injection of electrons into the first mini- 
band can be neglected. Therefore the current in the first 
mini-band vanishes near the boundary with the high-field 
domain and increases away from it, owning to the relax- 
ation from the second miniband. 

In doped superlattices the screening length is about 
one period. In other words, the drop of the field at the 
boundary between regions II and III causes the deple- 
tion of only the one well closest to the high-field domain 
and all other wells in the region III can be considered 
to be electroneutral. The relaxation of electrons from 
the second miniband, where their mobility is high, into 
the first miniband, where the electron mobility is much 
smaller, leads to the fields inhomogeneity on the scale of 
the relaxation length. We calculate the field distribution 
in the most interesting case, when this length is much 
larger than the superlattice period. 

A relaxation length large compared to the superlattice 
period and to the screening length allows us to use the 
condition of electroneutrality and to replace difference 
equations with differential ones. The total current j is 
the sum of the currents in the first mini-band, ji , and in 
the second mini-band, j2 '■ 

(5.26a) 
(5.26b) 



] = ]l+]2 , 

'^i f (1) TP rr.dn^^'' 



'-'2 I (2) TP rr 



dy 



(5.26c) 
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where fii and fj,2 are the mobihties in the first and the 
second miniband respectively. Here we assume the Boltz- 
mann distribution in both mini-bands and neglect the 
field dependence of the mobilities. The first assumption 
is reasonable because one can expect a significant heating 
of electrons injected into region III. The second assump- 
tion is natural due to the following circumstance. The 
field dependence of the mobility is appreciable when the 
field is close enough to the instability threshold field Fth. 
Usually the field in the high field domain Fh is signif- 
icantly larger than i^th- Equation ( S.lQj ) shows that in 
this case j* is significantly smaller than jo ^nd therefore 
Foo cannot be close to Fth- In the part of region III where 
the current goes mostly in the second miniband, the field 
is even smaller. 

The next two equations are the conditions of the elec- 
troneutrality and the continuity equations are 



dx 



n 



(2) 



T2ll 



(5.26d) 
(5.26e) 



In t he gen era l case the electric field can be found from 
Eqs. (|5.26b| ) - ( ^.26dD 



F 



J 



M2 J 



1-^ 



J2 

j 



(5.27) 



where Foo is the electric field far from the domain. Af- 
ter eliminat ing th e current j2 and the concentration n*-^^ 
from Eqs. ( 5.26c ), ( 5.26e ), and ( 5.27 ) one can get the 
differential equation for the field in region III 



r2 



F d F 



F 

Fqo 



1 



(5.28) 



The redistributing of electrons between two minibands 
and the field profile are characterized by two lengths: a 
diffusion length L-p and a drift length L^^, where 



7^2 



f^2T2lT 



and 



Ldt — f^2T2lFa 



(5.29) 



(5.30) 



The former length is the distance that electrons diffuse 
in the second miniband before the relaxation to the first 
miniband. The latter length is the distance that electrons 
run in the second miniband under the electric field be- 
fore the relaxation. Equation ( ^.28 ) has to be solved near 
the boundary for an arbitrary relation between these two 
lengths. Far from the boundary we calculate the distri- 
bution of electric field separately in two limiting cases, 
when one of these two lengths is much larger than the 
other. 

At y — there is no current in the first miniband and 



Eq. (5.27) gives the boundary condition for Eq. (5.2 



F 



x=0 M2 



(5.31) 



One can see that near the boundary F/Foa is small if 
Ml ^ M2. That can be expected because the second mini- 
band is usually wider than the first miniband. Thus we 
can neglect F/Fao on the right hand side of Eq. (5.28). 
This simplification allows us to solve this equation in 
terms of Airy functions. The result is 



F 

'fZ. 



r2 

7"2 



1/3 



Ai^ {{y - ya)lLa) 
Ai((y-2;o)/io) 



(5.32) 



where Ai'(^) denotes the derivative of the Airy function 
with respect to its argument and Lq — {2L'^/ Lj^-cY/^ . 

The other parameter yo can be found by substitution of 
Eq. ( 5.32| ) into the boundary condition Eq. ( |5.3l| ). This 

gives for ^ = -yo/Lo 



q - -Ai'(0/Ai(0 , 



(5.33) 



where q = (Mi/M2)(^nr/4i^^ 



,1/3 



The solution Eq. ( |5.32| ) shows that the electric field in- 
creases away from the boundary. Far from the boundary 
where F ^ Foo the solution Eq. (5.32) is not valid any 
more. In the case of a short drift length LDr^Cip o ne can 
neglect the first term on the right hand side of Eq. (5.2S) 



because it contains a small parameter Ldi/^c- We can 
show this by replacing y / L-p by a dim ensionless variable. 
In this case the solution to Eq. (5.28) is 
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^ oo 



Jl 
j 



Ml J2 
M2 j 



= 1 - 



1 - 



Ml 

M2 



cxp{-y/Lv) , 



(5.34) 



where the point a; = corresponds to the boundary be- 
tween regions II and III. The characteristic length of the 
electron relaxation from the second miniband to the first 
one is in this case the diffusion length L-p. 



In the opposite case when Lur^L-p, Eq. ( 5.28[ ) can 
be simplified in the region where [Ld / Lj^^^/'^ -^F / Fqq, 
namely, the term on the left hand side of Eq. ( 5.28 ) can 
be neglected. The solution to the resulting equation is 



F/Foo - log(l - F/Foo) ^{y- yo)/Lur 



(5.35) 



Here yo is the same number as in Eq. (5.32). This can 
be proven by matching the asymptotes of the both so- 
lutions, Eq. ( ^.32D and Eq. ( |5.35D , in the region where 

{Lv/LBr-P^<^F/Foo « 1. 

The very important property of Eqs. (5.26) is that the 
number of electrons in the second miniband explicitly 
depends on t he int erminiband relaxati on tim e. On e can 
get from Eq. ( ^.32D together with Eqs. ( ^.26e| ) and ( |5.27| ) 



(5.36) 



The parameter q, which defines the number of electrons 
in the second mini-band near the boundary with region 
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FIG. 5.5. The concentration of the electrons in the sec- 
ond mini-band in region III of the superlattice close to the 
boundary with region II is plotted versus the parameter 



1/3 

II, is proportional to . When q is small ^ « — 1 and 
7i(^^(0) w 2qn. When q is large the asymptote of the 
Airy function gives q « ^~"'^/4 -|- ^""^^^ and n2 ~ n. The 
concentration n^^^ as a function of q is plotted in Fig. 5.f; . 



The obtained field distribution is not valid at the 
distance of about one screening length from the high- 
field domain because the electro-neutrality condition 
Eq. (5.26d) cannot be used there. Thus in the well at 
the boundary between regions II and III, see Fig. 5 



where the field changes significantly, the concentration 
of electrons has to be found from the Poisson equation. 
The field change near the next well, i> = i>ii + 1 is 
much smaller and we can consider it and the rest of 
region III as electro-neutral. Therefore the solution to 
Eqs. (5.26) is valid up to well number = vu + 1, and 

"lii+i ^ n^'^^\x=o- The field in the barrier between wells 
number vu and vn -\- 1 can be easily computed from 
Eqs. ( ^.4a ) and (5.8t) with v = vn. In the case when 
the relaxation length of electrons from the second mini- 
band to the first one is much larger than one superlattice 
period we can neglec t the current between first levels and 
also put in the Eq. ( 5.4a ) jn =0 and ^22 = j- 



5.7. DISCUSSION 

In this section we give a qualitative description of 
the current-voltage characteristic of the entire superlat- 
tice. Under a small external bias the superlattice ex- 
hibits Ohm's law. At a higher bias the current reaches 



a maximum value. The uniform potential distribution 
in the superlattice is unstable at this point of the I-V 
characteristic. The thea[ai-|Of instability was developed 
by the present authorsE3Ej. The instability eventually 
leads to the formation of t he h igh-field domain and a 
sharp current drop. In Sec. 5.2 we argued that due to 
the charge accumulation at the high-field domain bound- 
ary, its length can be more than one superlattice period. 
Now we estimate this length. 

The threshold value of the bias just before the insta- 
bility point is NFtYil, where A'^ is the number of the su- 
perlattice periods. Right after the instability point this 
bias is a sum of the voltage drop across the high field 
region FuNul and the low field region w (A — Nu)Fool- 
Here we neglect a field inhomogeneity near the domain 
boundaries. So we have 



AFthZ = FhAiiZ + (A - Aii)FooZ , 



and 



A^ii = A^ 



Fth — Fa 

Fyi — Fry. 



(5.37) 



(5.38) 



This equation shows that in very long superlattices the 
length of the high-field domain is proportional to the 
length of the superlattice. This fact is a direct conse- 
quence of the upstream diffusion current near the bound- 
ary of the high-field domain. Without the diffusion the 
minimal length of the high-field domain is one superlat- 
tice period and Fth — ^00 c< 1/A^. Due to the diffusion 
current F^o is limited from above by a value in depe ndent 
of the superlattice length. From Eqs. ( ^.9a| ), ( 5.18| ), and 
( ^.19 ) we have 



F„ > Fa 



(5.39) 



For Fth > foe, one can easily see that Nn decre ases m ono- 
tonically with increasing Foo. Then Eqs. ( ^.38D and 



(^.39 ) give the following condition for A^n 



A^ii > A 



i^H + 3Fth - 



3(Fh- 



"thJ 



(5.40) 



Actually the field in regions I, II, and III is inhomoge- 
neous. As a result there is a correction to the right-hand 
side of Eq. ( ^.40 ). Usually this correction can be ne- 
glected because it does not contain the factor A^. 

With a further increase of the bias the current nearly 
periodically increases and drops down. After Esaki and 
Chang,El, each oscillation is associated with the extension 
of the high-field domain by one period and the number of 
oscill ation A"osc is expected to be equal to A^ — 1 . Equa- 
tion (5.40) shows that Aqsc can be less than A^ — 1 and 
that really is the|-Gase in some experiments. For instance, 
Kawamura et al.cJ observed a current- voltage character- 
istic with A^osc = 35 in a superlattice that had A^ = 39 
barriers, see Fig. 5.6. The difference between A^ and A'osc 
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FIG. 5.6. Well number dependence of oscillation repetition 
iber from the experimental work of the Y. Kawamura. et 



al 



can be interpreted as the formation of a high-field do- 
main with the minimal length Nn = 4. The value of the 
threshold field, i^th^ in this experiment can be obtained 
by dividing the threshold voltage, f^O.SV, by the number 
of supcrlattice barriers. The voltage drop across a barrier 
in the high- field domain F-al is equal to the period of the 
oscil lation s, 0.14V. The substitution of these values into 
Eq. Jslcl ) gives A'li > 3.1 which is in agreement with the 
value obtained from the number of the oscillations. 

In order to understand how a large domain comes 
about, we performed a numeric calculation of the insta- 
bility development. The simulation was made for 20 pe- 
riod superlattice with the inter-subband energy separa- 
tion lOOmeV, and the electron concentration correspond- 
ing to the Fermi energy 40meV. We assume some reason- 
able relations between transition and relaxation times. 
The parameters characterizing screening in the superlat- 
tice were \\ — 4, and \\ = 8. The results of the simu- 



lation are shown in Figs. 5.7 and 5.8. One can see that 
the domain starts to grow as a large scale instabilityt£l 
which transforms into a short domain with a length of 
three periods. 

In the case of low doping there are not enough electrons 
to form the depletion layer downstream of the doajiain. 
As the result the domain is formed near the anodeE3 and 
region III does not exist. The position of the high- field 
domain in a highly doped superlattice is determined by 
unintentional nonuniformity of the superlattice or by a 
fluctuation that initiated the domain. 

For fields larger than i^th , the oscillating current is lim- 
ited from above by the value of j* . It is important to note 
that this value is usually smaller than joj i- e. the current 



eF^d{t) 




FIG. 5.7. The time evolution of the voltage drops distribu- 
tion. The superlattice has 20 periods, cathode and anode are 
marked by C and A respectively. The interlevel spacing in the 
well is lOOmeV, whereas the domain height is only 70meV per 
period. 




FIG. 5.8. The time evolution of the concentration at the 
second level in each well. The superlattice has 20 periods, 
cathode and anode are marked by C and A respectively. One 
can see the tail of the electron distribution in region III (be- 
hind the domain). 
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value a,t F = Kv,, see Sec. 5.4. Such a feature of t 



V characteristic is typically observed in experimentcjH. 
We will show now that the high-field domain expands 
by joining another period from region I even in the case 
when region III does exist. For t his p urpose we make 
use of an equation similar to Eq. ( 5.37 ), where we take 
into account the nonuniformity of the field at the accu- 
mulation layer near the domain boundary. So instead of 
Eq. ( 5.37D we have 



= ((iVi - l)Foo -I- + iViiFff)? + l/m , (5.41) 

where V is the applied bias, N\ and A^n are the numbers 
of periods in regions I and II respectively, vi is the num- 
ber of the first barrier in the high-field doma in, V in is the 
potential drop across region III. Equation (5.41) resem- 
bles Eq. (5.37), however, it explicitly takes into account 
the field inhomogeneity in regions I and III. 

As the applied voltage increases, dV/dj goes to infin- 
ity which eventually leads to a current discontinuity and 
extension of the high field domain by one period. One 
can see from Eq. (5.41) that dV/dj — s- oo when in one 
of the barriers dF/dj oo. This cannot take place in 
region II, where the maximum value of the current is 
larger than that in region I. The infinite value of dF/dj 
or zero value of dj/dF is ultimately connected with the 
maximum in the I-V characteristic of one barrier. This 
maximum is reached only in the regions where current 
is confined in the first mini-band, i.e. in region I and 
in a part of region III. The largest field in these re- 
gions is in the las t barrier of region I, F^j^i, see the 



end of the Sec. 5.4. The field in this barrier correspond- 
ing to dF/dj — *■ oo is larger than Fth for the following 
reason. In this barrier dj/dF = djcond/dF + dj^is/dF. 
Since the field in region I grows faster than in region II, 
the concentration in the well between these regions de- 
creases. Therefore the diffusion current also decreases 
which means that djdm/dF > because the diffusion 
current is directed against the conduction current. As a 
result dj/dF = not when djcond/dF — (i.e. F = Fth), 
but when djcond/dF < 0, (i.e. F > Fth). 

The field in region II is about Eg /el and the field 
in region I is smaller or about Fth = F/eL Usually 
F is only a few pfr4- So, in experiments with a big 
energy separationcJ^EJ, Eg w 200meV, the field differ- 
ence between regions I and II is also large and it causes 
a large charge accumulation at the boundary between 
these regions. In these experiments eF^^l « 140 mcV , 
eFth^ « 12.5meV, so Fth/Fn ~ 0.089 and Eq. ^1^ ) 
gives j* ~ 0.30jo- That value is in good agreement with 
the ratio of the maximum of the current in the oscillating 
region to the peak value of the current at the instability 
threshold, see Fig. 2(b) in Ref. |2^, reproduced in Fig. 1.2 
of the present work. 

The oscillations of the current in a superlattice with a 
high-field domain caused by the expansion of this domain 



with an increase of the applied bias have a period which 
can be associated with the potential drop per superlattice 
period in the high-field region F^il. However, Eq. (5.24) 
shows that the potential drop per period in the domain 
is less than the energy separation between levels Eg by a 
significant quantity, which can reach a few tens of meV. 

Such a big difference between the period of the oscil- 
lations of the current-voltage characteristic and the en- 
ergy separation between the leiiels was detected in the 
experiment of Kawamura et al.til, see Fig. 3 in the cited 
work. One can see that the difference between the inter- 
subband energy space and the electric current oscillation 
period increases with Eg. This has a simple physical 
meaning: in samples with larger Eg the upper limit of 
the current j* is lower and therefore the resonance in re- 
gion II is weaker. Other examples of such a difference 
can be found in Refs. |8^, and 

The main assumption which was made in our calcula- 
tion concerns a small deviation of the electron distribu- 
tion function from the Fermi function with effective tem- 
perature and chemical potential. Necessary conditions 
justifying this assumption (see Sec. Ill) may not be sat- 
isfied in the high-field domain. However, the features of 
the current-voltage characteristic discussed in the present 
work do not depend on the detailed shape of the electron 
distribution. 

In conclusion, we have studied the field distribution 
and I-V characteristic of the superlattice in the voltage 
region when a high- field domain exists. The accumula- 
tion of electrons at the domain boundary causes a strong 
limitation of the current. Due to this limitation the min- 
imal length of the high-field domain can be not one but 
a few superlattice periods. The current limitation also 
results in the reduction of the period of the current os- 
cillations compared to that corresponding to the energy 
separation between the first and second levels in a well. 
These results are in good agreement with available exper- 
imental data. The field distribution in the region down- 
stream of the high-field domain (if it exists) is nonuniform 
due to electron relaxation from the second miniband to 
the first miniband. 



References 



^ cited after Michel De Montaigne, Essays (Penguin books, 
Baltimore, 1967). 

^ C. Wcisbuch, in Application of Multiquantum Wells, Selec- 
tive Doping and Superlattices., Vol. 24 of Semiconductors 
and Semimetals, edited by R. Dingle (Academic Press Inc, 
San Diego, 1987), Chap. I, p. 1. 

^ L. Esaki and R. Tsu, IBM J. Res. Dev. 14, 61 (1970). 

* R. F. Kazarinov and R. A. Suris, Fiz. Tekn. Poluprovodn. 

5, 797 (1971), [Sov. Phys.-Semicond. 5, 707 (1971)]. 

^ R. F. Kazarinov and R. A. Suris, Fiz. Tekn. Poluprovodn. 

6, 120 (1972), [Sov. Phys.-Scmicond., 6, (1972)]. 

'5 L. Esaki and L. L. Chang, Phys. Rev. Lett. 33, 495 (1974). 

F. Capasso, K. Mohammed, A. Y. Cho, IEEE J. Quantum. 
Electron. QE-22, 1853 (1986). 

* B. F. Levine, K. K. Choi, C. G. Bethea, J. Walker, and 
R. J. Malik, Appl. Phys. Lett. 51, 934 (1987). 

" M. Helm, P. England, E. Colas, F. DeRosa, and S. J. Allen, 

Jr., Phys. Rev. Lett. 63, 74 (1989). 
^° J. Faist, F. Capasso, C. Sirtori, D. L. Sivco, J. N. Bail- 

largcon, A. L. Hutchinson, S.-N. G. Chu, and A. Y. Cho, 

Appl. Phys. Lett. 68, 3680 (1996). 

G. Schwarz, A. Wacker, F. Prengel, E. Scholl, J. Kastrup, 

H. T. Grahn, and K. Ploog, Semicond. Sci. Technol. 11, 
475 (1996). 

H. Le Person, C. Minot, L. Boni J. F. Palmier, and F. Mol- 
lot, Appl. Phys. Lett. 60, 2397 (1992). 
"S. H. Kwok, T. B. Norris, L. L. Bonilla, J. Galan, 

F. C. Martinez, J. M. Molera, H. T. Grahn, K. Ploog, and 
R. Merlin, Phys. Rev. B 51, 10171 (1995). 

H. Grahn, J. Kastrup, K. Ploog, L. Bonilla, J. Galan, 
M. Kindelan and M. Moscoso, Jpn. J. Appl. Phys. 34, 4526 
(1995). 

1^ R. Klahn, S. H. Kwok, H. T. Grahn, and R. Hey, Phys. 

Rev. B 52, 8680 (1995). 
^'^ A. Wacker, M. Moscoso, M. Kindelan, and L. L. Bonilla, 

Phys. Rev. B 55, 2466 (1997). 
" J. Kastrup, R. Hey, K. H. Ploog, H. T. Grahn, L. L. BoniUa, 

M. Kindelan, M. Moscoso, A. Wacker and J. Galan, Phys. 

Rev. B 55, 2476 (1997). 

A. Sibille, J. F. Palmier, and F. Mollot, Appl. Phys. Lett. 
60, 457 (1992). 

K. Fujiwara, K. Kawashima, T. Yamamoto, N. Sano, 
R. Cingolani, H. T. Grahn, and K. Ploog, Phys. Rev. B 
49, 1809 (1994). 

A. Sibille, in Semiconductor Superlattices, Growth and 
Electronic Properties, edited by H. T. Grahn (World Sci- 
entific, London, 1995), Chap. 2. Miniband transport. 

G. Bastard, J. A. Brum, and R. Ferreira, in Semiconductor 
Heterostructures and Nanostructures, Vol. 44 of Solid State 
Physics, edited by H. Ehrenreich and D. Turnbull (Aca- 



demic Press, Inc, San Diego, 1991), Chap. II. The envelop 
function approximation, p. 232. 

22 D. L. Smith and C. Mailhiot, Rev. Mod. Phys. 62, 173 
(1990). 

23 B. Laikhtman, Phys. Rev. B 44, 11260 (1991). 

2* Y. Kawamura, K. Wakita, H. Asahi, and K. Kurumada, 
Jap. J. Appl. Phys. 25, L928 (1986). 

2^ Semiconductor Superlattices, Growth and Electronic Prop- 
erties, edited by H. T. Grahn (World Scientific, London, 
1995). 

2*5 B. Laikhtman and D. Miller, Phys. Rev. B 48, 5395 (1993). 
2'^ D. L. Miller and B. Laikhtman, Phys. Rev. B 50, 18426 

(1994) . 

28 S. H. Kwok, R. Merlin, H. T. Grahn and K. Ploog, Phys. 

Rev. B 50, 2007 (1994). 
2^^ J. Kastrup, H. T. Grahn, K. Ploog, F. Prengel, A. Wacker, 

and E. Scholl, Appl. Phys. Lett. 65, 1808 (1994). 
3° F. Aristone, J. F. Palmier, P. Gassot, J. C. Portal, and 

F. Mollot, Appl. Phys. Lett. 67, 2916 (1995). 

31 E. M. Epshtein, Izv. Vuzov, Radiofizika 22, 373 (1979). 

32 F. G. Bass, A. A. Bulgakov, and A. P. Tetervov, High- 
frequency properties of semiconductors with superlattice 
(Novoscience Publisher Inc., New- York, in press 1996), 
Chap. 7. Conduction and galvanomagnetic effects. 

33 D. L. Miller and B. Laikhtman, Phys. Rev. B 52, 12191 

(1995) . 

3* I. B. Levinson and Ya. Yasevichute, Sov. Phys.- JETP 35, 
991 (1972). 

3^ R. A. Suris and B. S. Shchamkhalova, Sov. Phys. Semicond. 

18, 738 (1984). 

36 M. Lee, N. S. Wingrcen, S. A. Solin, and P. A. Wolff, Solid 
State Comm. 89, 687 (1994). 

3^ D. L. Miller and B. Laikhtman, Phys. Rev. B 54, (1996). 

3* E. E. Mendez, F. AguUo-Rueda, and J. M. Hong, Phys. 
Rev. Lett. 60, 2426 (1988). 

3^ K. Kawashima, T. Yamamoto, K. Kobayashi, and K. Fuji- 
wara, Phys. Rev. B 47, 9921 (1993). 

G. Cohen, I. Bar-Joseph, and H. Shtrikman, Phys. Rev. B 
50, 17316 (1994). 

*^ R.Tsu and G. Dohler, Phys. Rev. B 12, 680 (1975). 

''2 P. Leisching, P. Haring Bolivar, W. Beck, Y. Dhaibi, 

F. Briiggemann, R. Schwedler, H. Kurz, K. Leo, and 

K. Kohler, Phys. Rev. B 50, 14389 (1994). 
*3 G. Von Plessen, T. Meier, J. Feldmann, E. O. Gobel, 

P. Thomas, K. W. Goosen, J. M. Kuo, and R. F. Kopf, 

Phys. Rev. B 49, 14058 (1994). 

R. Peierls, Z. Physik 80, 763 (1933). 

J. M. Luttinger, Phys. Rev. 84, 814 (1951). 
'^^ E. I. Blount, Phys. Rev. 126, 1636 (1962). 

A. M. Berezhkovskii and R. A. Suris, Zh. Eksp. Teor. Fiz. 



53 



5.7. DISCUSSION 



54 



86, 193 (1984), [Sov. Phys. ,]ETP 59, 109 (1984)]. 

A. A. Abrikosov, Fundamentals of the theory of metals 
(North-Holland, Oxford, 1988), Chap. 10 p. 198. 

"'^ E. M. Lifshitz and L. P. Pitaevsky, Statistical Physics, 
Part 2, Vol. 9 of Landau and Lifshitz: course of theoret- 
ical physics (Pergamon Press, Oxford, 1980), Chap. VI. 
Electrons in the crystal lattice. §56. 

L. Esaki, in Molecular Beam Epitaxy and Heterostructures, 
edited by L. L. Chang and K. Ploog (Martinus Nijhoff Pub- 
lishers, Dordrecht, 1985), Chap. 1. 

51 K. K. Choi, B. F. Levine, N. Jarosik, J. Walker, and 
R. J. Malik, Phys. Rev. B 38, 12362 (1988). 

52 T. H. H. Vuong, D. C. Tsui, and W. T. Tsang, J. Appl. 
Phys. 66, 3688 (1989). 

53 R. A. Davies, D. J. Newson, T. G. Powell, M. J. Kelly, and 
H. W. Myron, Scmicond. Sci. Technol. 2, 61 (1987). 

5* A. Sibille, J. F. Palmier, A. Celeste, J. C. Portal, and 

F. MoUot, Europhys. Lett. 13, 279 (1990). 

55 M. Lee, S. A. Solin, and D. R. Hines, Phys. Rev. B 48, 
11921 (1993). 

5*^ W. Miillcr, H. T. Grahn, K. von Klitzing, and K. Ploog, 

Phys. Rev. B 48, 11176 (1993). 
5^F. Aristone, A. Sibille, J. F. Palmier, D. K. Maude, 

J. C. Portal, and F. MoUot, Physica B 184, 246 (1993). 
5* H. J. Hutchinson, A. W. Higgs, D. C. Herbert, and 

G. W. Smith, J. Appl. Phys. 75, 320 (1994). 

5^ A. Ya. Shik, Fiz. Tekh. Poluprovodn. 7, 261 (1973), [Sov. 

Phys.- Semicond. 7, 187 (1973)]. 
^° T. Ando, J. Phys. Soc. Jpn. 50, 2978 (1981). 

R. A. Suris and B. S. Shchamkhalova, Fiz. Tekh. 

Poluprovodn. 24, 1638 (1990), [Sov. Phys. Semicond. 24, 

1023 (1990)]. 

^2 B. Movaghar, Semicond. Sci. Technol. 2, 185 (1987). 

,J. F. Palmier, A. Sibille, G. Etcmadi, A. Celeste, and 
J. C. Portal, Semicond. Sci. Technol. 7, 283 (1992). 

^ R. Z. Sagdeev, D. A. Usikov, and G. M. Zaslavsky, Non- 
linear Physics (Harwood Academic Publishers, New York, 
1988). 

•^5 E. M. Lifshitz and L. P. Pitaevsky, Physical kinetic, Vol. 10 
of Landau and Lifshitz: course of theoretical physics (Perg- 
amon Press, Oxford, 1980), Chap. HL Landau damping 
§30. 

T. Ruf, J. Spitzcr, V. F. Sapcga, V. I. Belitsky, M. Cardona, 
and K. Ploog, Phys. Rev. B 50, 1792 (1994). 
^"^ L. L. Chang, H. Sakaki, C. A. Chang, and L. Esaki, Phys. 
Rev. Lett. 38, 1489 (1977). 

E. M. Lifshitz and L. P. Pitaevsky, Physical Kinetics, 
Vol. 10 of Landau and Lifshitz: course of theoretical physics 
(Pergamon Press, Oxford, 1980), Chap. IX. Metals. §84. 
Yu. M. Galperin, V. D. Kagan, and V. I. Kozub, Zh. Eksp. 
Teor. Fiz. 62, 1521 (1972), [Sov. Phys.- JETP 35, 798 
(1972)]. 

™ I. S. Gradsteyn and I. M. Ryzhik, Tables of Integrals, Series 
and Products (Academic, New York, 1980). 

B. Laikhtman and E. L. Altshuler, Ann. Phys. 232, 332 
(1994). 

P. G. Harper, Proc. Phys. Soc. London, Sect. A 68, 879 
(1955). 

^3 L. B. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1525 (1964), [Sov. 
Phys. JETP 20, 1018 (1965)]. 



^* E. M. Lifshitz and L. P. Pitaevsky, Physical Kinetics, 
Vol. 10 of Landau and Lifshitz: course of theoretical physics 
(Pergamon Press, Oxford, 1980), Chap. X. §93, §94. 

■^5 S. V. Gantsevich and V. L. Gurevich, Fiz. Tverd. Tela 7, 
2400 (1965), [Sov. Phys. Solid State 7, 1939(1966)]. 

'■'^ V. L. Gurevich and R. Katilius, Zh. Eksp. Teor. Fiz. 49, 
1145 (1965), [Sov. Phys. JETP 22, 796(1966)]. 
B. J. Skrommc, R. Bhat, M. A. Koza, S. A. Scliwarz, 
T. S. Ravi, and D. M. Hwang, Phys. Rev. Lett. 65, 2050 
(1990). 

''^ F. Piazza, L. Pavesi, H. Cruz, M. Micovic, and C. Mendoga, 

Phys. Rev. B 47, 4644 (1993). 
''^ R. Landauer, in Electrical Transport and Optical Properties 

of Inhomogeneous Media, edited by J. C. Garland and D. B. 

Tanner (AIP, New York, 1978). 
*° D. Bimberg, J. Christen, T. Fukunada, H. Nakashima, 

D. E. Mars and J. N. Miller, J. Vac. Sci. Technol. B. 5, 

1191 (1987). 

*i Y. Kawamura, K. Wakita, and K. Oe, Jap. J. Appl. Phys. 

26, L1603 (1987). 
82 K. K. Choi, B. F. Lcviuc, R. J. Malik, J. Walker, and 

G. G. Bethea, Phys. Rev. B 35, 4172 (1987). 
*3 B. F. Levine, K. K. Choi, C. G. Bethea, J. Walker, and 

R. J. Malik, Appl. Phys. Lett. 50, 1092 (1987). 

K. K. Choi, B. F. Levine, C. G. Bethea, J. Walker, and 

R. J. Malik, Appl. Phys. Lett. 50, 1814 (1987). 
85 H. T. Grahn, H.Schneider, and K. von Klitzing, Phys. Rev. 

B 41, 2890 (1990). 

M. Artaki and K. Hess, Phys. Rev. B 37, 2933 (1988). 
*^^D. C. Herbert, Scmicond. Sci. Technol. 3, 101 (1988). 
V. I. Falko, Phys. Rev. B 47, 13585 (1993). 
F. Prengel, A. Wacker, and E. Scholl, Phys. Rev. B 50, 
1705 (1994). 

^° C. H. Yang and S. A. Lyon, Physica 134B, 305 (1985). 
^1 A. Seilmeier, H.-J. Hiibner, G. Abstreiter, G. Weimann, 

and W. Schlapp, Phys. Rev. Lett. 59, 1345 (1987). 
^2 J. K. Jain and S. Das Sarma, Phys. Rev. Lett. 62, 2305 

(1989). 

^3 R. J. Bauerlc, T. Elacsscr, W. Kaiser, H. Lobentanzer, 
W. Stolz, and K. Ploog, Phys. Rev. B 38, 4307 (1988). 
M. C. Tatham, J. F. Ryan, and C. T. Foxon, Phys. Rev. 
Lett. 63, 1637 (1989). 

"5 H. T. Grahn, H. Schneider, W. W. Riililc, K. von Klitzing, 
and K. Ploog, Phys. Rev. Lett. 64, 2426 (1990). 
D. Y. OberU, D. R. Wake, M. V. Klein, J. Klem, T. Hen- 
derson, and H. Morkog, Phys. Rev. Lett. 59, 696 (1987). 
B. Deveaud, A. Chomcttc, B. Lambert, and A. Regreny, 
Solid State Gommun. 57, 885 (1986). 

^8 H. T. Grahn, Physica Scripta. T49, 507 (1993). 



